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(JZ i Abstract: We prove diffusion for a quantum particle coupled to a field of bosons (phonons or photons). The 
. _ ' importance of this result lies in the fact that our model is fully Hamiltonian and randomness enters only via the 
^vq . initial (thermal) state of the bosons. This model is closely related to the one considered in [5] but various restrictive 
assumptions of the latter have been eliminated. In particular, depending on the dispersion relation of the bosons, 
r 3 1 . the present result holds in dimension d > 3. 

-£ 1 Introduction 

c3 : 

The rigorous derivation of long-time diffusion from first principles of mechanics, be it quantum or classical, remains 
an inspiring challenge in mathematical physics. To our best knowledge, there are up to this date very few results of 
this type, see Section 12.51 for a brief review. Recently, in [5] , a model was introduced which is quite tractable and 
for which diffusion was proven in dimension d > 4. It is a quantum system described by a Hamiltonian of the type 

fTj ; H = H S + H E + \H I , Ael (1.1) 

oo ■ 

where ffs is the Hamiltonian of a free particle moving on the lattice and it consists of two parts ffs = ffkm + ffspin 
i^l describing the translational (kinetic) degrees of freedom and a spin degree of freedom, respectively. The Hamiltonian 

^^ ■ He describes a free field of bosons (the environment), and ffj effectuates the coupling between both. The system 
t— I , is started with the environment in a thermal state at inverse temperature j3 or in a nonequilibrium state (then we 

have two phonon fields, at different inverse temperatures f3\ ^ $2). Such models are a paradigm of open quantum 
J> , systems. The form of the Hamiltonian will be given in Section [21 but let us already list the properties that allow us 

to handle this model: 



• The mass of the particle, or, since we are on a lattice, rather the inverse hopping strength, is chosen large. 



• 



In (jl.ip this is accomplished by choosing -ffkin in ffs small. This allows a better control of a diagrammatic 
expansion in real space, since the particle needs a long time to explore a large volume on the lattice. 

Even though the mass is large, the 'mixing rate' that the spin and momentum degrees of freedom of the 
particle experience due to the interaction with the phonons is not small. This is possible because of the 
inclusion of the spin-degree of freedom. 
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• By choosing the interaction Hamiltonian sufficiently smooth (in the momentum of the phonons) and the 
dimension d sufficiently large, we ensure that the free space-time correlation functions of the boson field 
decay at an integrable rate in time. We need them to decay at least as 0(t~( 1+a )) for large times t with 
a > 1/4. To engineer this in d > 3, it suffices to choose the dispersion relation of the bosons to be quadratic 
in the momentum for small momenta, and to cut off the interaction for large momenta. 

• By choosing the coupling constant A small, we have a well controlled Markovian approximation (Lindblad 
equation) that describes the particle for times of 0(A~ 2 ). This Markovian approximation serves as a first 
approximation to the true behavior and we set up an expansion to control the deviations from it. 

Under these assumptions we prove that the reduced dynamics of the particle is diffusive. Our proof is based 
on a renormalization group (RG) method that was developed in [4j [5l H] to prove diffusion for random walk in 
a random environment (RWRE). In the present context the random environment is provided by the phonon field. 
Unlike in the case of RWRE, in the case at hand the particle influences the environment and the reduced dynamics is 
non-Markovian. However, the Markovian approximation mentioned above provides a starting point for the analysis 
where a Markovian dynamics is perturbed by a small non-Markovian noise. In units of the weak coupling time scale 
0(A~ 2 ) our model can then be viewed as a (quantum) random walk in a (quantum) random environment. The RG 
method consists of an iterative scheme to show that on successive larger temporal and spatial scales the random 
environment becomes smaller and smaller and the dynamics tends to a renormalized Markovian " fixed point" . We 
show that the renormalized noise vanishes in this limit by showing that its (quantum) correlation functions tend to 
zero. Here we use a formalism developed earlier by us [22] for the confined case, i.e. the proof that the state of a 
confined quantum system interacting with a similar field as here tends to the equilibrium state. 

The difference of the model considered in this paper and the one treated in [8] is that in the latter case an 
additional condition was imposed on the free boson correlation function that restricts the model to dimensions 
d > 4 and to a rather special class of analytic particle-phonon interaction terms. In the context of these models 
where the particle mass is chosen to scale as 0(A~ 2 ) it still remains a challenge to treat more generic phonon or 
photon reservoirs where the temporal correlations decay as 0(£ _1 ) (which is the case in d = 3 if the dispersion 
relation is linear for small momenta). To deal with these cases with our method one needs a more careful RG 
analysis. A much more difficult and interesting problem is to relax the large mass assumption. In this case the 
control of the corrections to the Markovian approximation seems still beyond current techniques. 

Acknowledgements. We thank for European Research Council and Academy of Finland for financial support. 

2 The model 

We consider a finite, discrete hypercube A = A^ = 7L d jUL d with LeN. We will take the thermodynamic limit at 
the end by letting L — > oo. Since most concepts in the present section depend on the volume A, we do not indicate 
it explicitly. 

2.1 Dynamics 

2.1.1 Particle 

The Hilbert space of the particle is 

,j#> s =y>(g>l 2 (A) (2.1) 

where S? is the finite dimensional space of internal degrees of freedom. The Hamiltonian is 

#S = # S pin <8> 1 + 1 <8> flidn, #k in = - — A (2.2) 

A z m v 

where A is the discrete Laplacian on A with periodic boundary conditions, defined by the kernel A(x, y) = 6\ x —y\ t i ~ 
2d8 x>v where \x — y\ — X)i=i \ x i ~ Ui\ with | • | the distance on the torus Z/XZ. The parameter m p is the mass of 



the particle, it is written with A~~ in front of it to remind us that we choose the mass very big by making A small. 
For simplicity, we assume that the Hamiltonian H spm is nondegenerate such that we can label its eigenvectors by 
their eigenvalues, which we denote by e £ cr(H sp i n ). 

2.1.2 Phonon field 

A single phonon is described by the one- particle Hilbert space ^ 2 (A). It will however be convenient to consider this 
space in the Fourier- representation, Z 2 (A*) where A* = ^-(Z d /LZ d ) is the dual lattice. We consider A* as a subset 
of T d -the d-dimensional torus identified with [—ir, 7r] d . The Hilbert space of the phonon field is 

^ E = r(Z 2 (A*)) (2.3) 

where r(f)) is the symmetric (bosonic) Fock space built on the one-particle Hilbert space f). 
The Hamiltonian of the phonon field is given by 

H E =J2 "(lK q a q (2.4) 

<?eA* 

where a*/a q are the creation/annihilation operators satisfying the canonical commutation relations (CCR) 

[aq,a*q>] = Kq'- ( 2 - 5 ) 

w (<?) > is the frequency of the phonon with momentum q and we take to a smooth function defined on T d . We 
impose later further conditions on u. 

We will also consider a non-equilibrium setup, in which case we consider two different phonon fields, distinguished 
by the label j = 1, 2. In this case the Hilbert space is J?k = T(/ 2 (A*)) (g) T(l 2 (A*)) and the Hamiltonian 

h e=^2 {^i{q)a* q ,ia q s + oj 2 {q)a* :2 a q ,2) ( 2 - 6 ) 

where a* i/a q ^i act on the first tensor in Jfy and a* 2 / a q ,2 on the second. 

We assume the reader is familiar with these basic notions of second quantization and we refer to 9 for more 
detailed accounts. 

2.1.3 Full Hamiltonian 

The interaction between particle and phonon field is chosen linear in the creation/annihilation operators. It is of 
the form 

Hi= J2 (W ® e iqX ®<P 3 (q)a qiJ ) +h.c. (2.7) 

geA*,i=l,2 

where 'h.c' stands for 'hermitian conjugate', W is a Hermitian matrix acting on ,5f (the space of the internal degree 
of freedom), X is the multiplication operator with the variable x £ Z d on l 2 (Z d ), and <f>j(q) are "structure factors" 
that describe the coupling to the phonons. We will take 4>j(q) the values of smooth functions <pj : T d — > C at q £ A*. 
We will impose some further conditions on <pj later. 

The fact that the particle couples to both phonon fields via the same W^-matrix, is by no means important, and 
we make it for notational simplicity. 

The total Hamiltonian is then 

H = H S + H E + Aifj (2.8) 

acting on J$? = Jtfc, Cg> J^,. If ujj(q) > for any q £ A*, then Hi is an infinitesimal perturbation (in the sense of 
Banach operator theory) of H E . By a standard application of the Kato-Rellich theorem, H is self-adjoint on the 
domain of He (we need not worry about H§ since it is bounded). The condition U)j(q) > is also necessary to have 
a well-defined finite-volume Gibbs state (see e.g. (|2.16jl ). However, once we take the thermodynamic limit, the only 
remaining infra-red regularity assumption will be Assumption lAl (in particular U)j(0) may vanish). 



2.2 States 

States of the combined system and environment are given by density operators pse G 38\(J%s J^e) where 3§\ 
denotes trace class operators. The initial states that we will consider are of the following form 

Pse = Ps ® pt (2.9) 

with 



tf* = l- c -^ H ^i ® e -P2H E2 



AA = TV (e _ftHEl 8e " ft ^) (2.10) 

and ps G 3B\{3tf§) is chosen to have support concentrated around the origin in the following sense. Elements of J^s 
can be represented by functions ip{x) with x E A C Z d and taking values in the spin space S^. In this basis ps is 
given by a kernel (matrix) p${x' ,x) taking values in 38(5^). Then we require 

p s {x',x)=0, for \x\,\x'\ >R (2.11) 

for some R. 

Since the environment Hamiltonians are quadratic in the creation/annihilation operators, the density matrix 
Pg is a "Gaussian state" , sometimes also referred to as a 'quasifree state" . Moreover, since these density matrices 
are functions of the environment Hamiltonian, the initial environment state p E is obviously invariant under the 
free environment evolution: 

e -itH Ep ref e itH E = ^ref ^^ 

In fact, these two properties, Gaussianity and invariance under the free dynamics, are what we will really use in 
our analysis, and the specific choice that we made for p^ is not important, except for one of our results which 
will additionally require /3i = 02 and where we exploit the fact that the system is in (close to) thermodynamic 
equilibrium. The dynamics of the density matrix of the entire system is given by 

PSB,t - c- itH (p s p^ f ) c itH = c- itL (p s ® p% { ) (2.13) 

where we denote L := ad(H). 

The reduced dynamics of the particle is defined by taking a partial trace Tr E over the environment Hilbert space 

JTe. 

PS,* - Tr E pse,* = Tr E (e- itL ( / 9 S p r E ef )) := Z tPs (2.14) 

All our results will concern the reduced density matrix ps,t, which means that we only consider particle observables, 
but it is straightforward to extend the formalism so as to handle observables of the boson field, as well as more 
general initial states. 

2.3 Thermodynamic limit 

Since we are ultimately interested in long-time properties, we have to perform the thermodynamic limit to eliminate 
Poincare recurrences. We take care of this below, but first we introduce the free boson correlation function which 
will play an important role in our our main assumptions. First, we define the so-called Segal field operators 

$(*,*):= Y, (e l{qX+ " j(q)t) $Mah+z~' l{qx+wMt) JMa<i,^ ( 2 -15) 

<reA*,j={i,2} 

The correlation function is then defined as 

COM) := Tr[p£ f $(.T,i)$(0,0)] 

= -L Y U.( g )|2f^_l e i(*«+to,(«)) + L_^-i(*«+**,(«)A (2.16) 

1 ' qeh',j={l,2} V 7 



where the second equality is again a standard exercise in second quantization, and one recognizes the Bose-Einstein 
distribution l/(e^ Wj — 1). The correlation function C,{x, t) will naturally come up in the evaluation of the perturba- 
tion series for the dynamics. 

To describe the thermodynamic limit, let us indicate the dependence in (|2.16[) on the volume A explicitly by a 
superscript ( (x,t). We then assume that the finite- volume free correlation functions converge 

lim C A (M) = C(M) (2.17) 

K/-7L d 

uniformly on compacts (in t £ K), and that the limiting correlation function £(2^) is bounded. This will hold for 
the basic examples of optical and acoustic phonons in Section 12.41 (in the latter case we need to exclude q = in 
the sum in eq. (|2.16|) . cfr. the discussion in Section r2.1.3|) . Note also that the spaces Jrf?s for A finite are naturally 
embedded into the Ji?s for A = Z d (L = oo). In particular, a density matrix pg satisfying (|2.11|) . is embedded into 
•^i(^s) for L > R. More generally, we have 

Lemma 2.1. If the initial state is chosen as described above (including in particular the convergence (|2.17|) and the 

boundedness of Q{x,t)) then the limit 

ps,t := Urn p£, 4 (2.18) 

Ayz d 

exists in 3§i(J#s) with J%s defined with A = 1/ . 

2.4 Results 

We need an assumption that will guarantee sufficiently fast decay of correlations of the freely evolving phonon field. 
Assumption A (Decay(a)). There is an a > such that 



dt (1 + \t\) a sup \C(x, t)\ < oo (2.19) 

xei, d 



With no loss we suppose a < 1. 



The next assumption eliminates a drift by requiring that the model be reflection-symmetric and rotation- 
symmetric (as far as the lattice allows). Let O^d be the subgroup of the orthogonal group O(d) that fixes the lattice, 
i.e. O G O zd iff. ueZ d ^>OueZ d and \Ou\ 2 = \u\ 2 with \u\l = J2i=i N 2 - 

Assumption B (Symmetries). The correlation function £ is 0%d -invariant in the sense that 

C(Ox,t) = ((x,t), for any O £ Z d. (2.20) 

Of course, this assumption is satisfied by choosing the parameters of the Hamiltonians i?E, Hi and the initial state 
;0g f symmetric. The Z d-invariance of the particle dynamics was already assured by our explicit choice of -ffkin- 

The next assumption assures that the particle is sufficiently well-coupled to the phonon field. To state it, we 
need some definitions. 

First, we assume that the spectrum of -ff sp in is non-degenerate. Then we can label a basis ip e € 5? of eigenvectors 
of -ff S pin by e £ cr(-ffspin) and we define W e ^i := (ip e , Wipe')j? (in general, we write (•, -)j? for the scalar product in 
a Hilbert space J^). 

Secondly, let us introduce the set of Bohr frequencies corresponding to the spin Hamiltonian if spin, 

£ := <r{&d{H spin )) = { £ = e-e'- e, e' G <r(H spin )} (2.21) 

We assume the spectrum of ad(-ff sp i n ) is also non-degenerate except for the eigenvalue which is dimJ^-fold 
degenerate. That is, for any e € £, e ^ 0, there is a unique pair (e, e') c cr(7J sp i n ) such that e' — e = e. Next, let 



Ux)= / dte- w X(x,t) (2.22) 



be the Fourier transform of the correlation function in time (well-defined by Assumption |A|) . We require that 

lim Ce(x) = (2.23) 

\x\— >oo 

for all ^ e E £. 

From (the infinite volume limit of) (I2.16P the distributional Fourier transform in x of £ w is formally 

Uq) = (-l)^M 27 r J2 \<t>M?-p^—cK"M |u,|). (2.24) 

.7=1,2 

We assume that for all ^ e E £ the distribution £ e defines a finite positive Borel measure on T d which we denote by 
Ce(dg). It is readily seen that this assumption holds e.g. if we suppose that V(jJj nowhere vanishes on a neighborhood 
of the set Ai e .j = {q E T d \ ujj(q) — \e\} for ^ e E £. For to — 0, we assume (^ = 0. These conditions hold in the 
examples below. 

We will now define a Markov jump process with state space <r{H sp i n ) x T d . The rates j(e',dk'; e, k) of the 
jump process (jumps are from unprimed to primed variables) are translation invariant in fc-space: j(e', dfc'; e, k) = 
j(e', d(k' — k); e, 0) and given by 

j(e', dfc'; e, 0) = \W e y \ 2 Q^ e (dk') (2.25) 

The main assumption then is that this Markov process is irreducible. We summarize: 

Assumption C (Fermi Golden Rule). Assume 

• The spectra of H sp i n and &d(H sp i n ) (apart from the eigenvalue 0) are non-degenerate 

• £o = and for all =/= e E <j(&d(H sp i n )), ( e (x) decays at infinity and £ e defines a finite positive measure. 

• The Markov process on o~(H sp i n ) x T d with jump rates \2. 25\) is irreducible. 

Before proceeding to the results, let us give some examples of models for which all our assumptions are satisfied. 
Example 1 (optical phonons). Let 

d 

cu(q) = (m 2 ph + ^ sin 2 q, ) i , with m ph ^ . (2.26) 

4=1 

This is a typical dispersion relation of an optical phonon branch. The form factor is (j>(q) is a smooth function on 
M. d which is nonzero only in neighborhood of such that |detHess(u;)| is bounded below on Supp</>. Then a standard 
argument relying on stationary phase estimates yields 

sup|C(^*)| < C(l + |t|)- d / 2 (2.27) 

X 

Hence Assumption [S] is satisfied with a = (d/2) — 1 — 5, for any 5 > 0. To ensure that Assumption [Cl is satisfied, we 
choose ,¥ — C 2 with W — a x and -ff sp i n = £0^2 where we use the traditional notation a x = ( -, ) , a z = I „ , J 
Then a sufficient condition for Assumption [Cl is that A4 eo HSupp0 has positive measure on Ai eo . With these choices, 
our equilibrium result (/3i = /?2) holds for d > 3, and the non-equilibrium result (/3i ^ P2) holds for d > 4. 

Example 2 (acoustical phonons). If one considers acoustical phonons, for example with the above dispersion relation 
(I2.26P with m p h = 0, then, for smooth </>(g), sup^, |C(a:,i)| < C(l + |t|) _ ^ d_1 ^ 2 . Hence, the equilibrium result holds 
for d > 4 and the nonequilibrium one for d > 5. From the point of view of our techniques, the important difference 
with the above example lies not therein that inf uj = 0, but in the fact that w(g) is linear in |g| for small q. 



2.4.1 Diffusion 

Our most important result concerns diffusion of the particle. To state this somehow concisely, we note that the 
density matrix ps,t determines a probability measure F t (x) on Z d : 

W t (x) = Tts'\ps,t(x,x)]= ^2 Ps,t(x,e;x,e) 

eecr(ff spin ) 

where Try is the partial trace on 33(5^) (below Tr is the full trace on Jf?s), and we wrote the matrix ps,t(x',x) G 
B{y) in the basis indexed by a(H sp i n ) as ps,t(x' ,e';x, e). Physically speaking, Pt(x) is the probability to find the 
particle at time t on the lattice site x E Z d . Equivalently, for an observable F(X) 

J2^t(x)F(x)=Ti[ps,tF(X)}, Fel°°(Z d ) (2.28) 

X 

Hence, the function 7 1— > Tr[ps it e nX ], figuring in the theorems below, is the characteristic function of the probability 
measure F t . 

Theorem 1 (Equilibrium). Assume that Assumptions [3 [0 and Assumption\^\with a > 1/4, hold, and moreover, 
that j3\ — /?2- Then, there is a Xq > such that, for < |A| < Ao, the characteristic function of -7= converges to a 
Gaussian: for all k £ M. d 

hmTr[e ife ^p st ]=e-l fc l 2 ^, (2.29) 

t— >oo 

for some strictly positive diffusion constant D* . Moreover, also moments of -*= converge to moments of the Gaussian. 
That is, for any multi-index I 

lim Tr[(4)W] = Hd fe )Vl fc l 2D * . (2.30) 

Theorem 2 (Non-Ecjuilibrium). Assume that Assumptions\B^\^ and AssumvtionVfi\ hold with a > 1/2, but possibly 
Pi 7^ h- Then, there is a Xq > such that, for < |A| < Ao, the conclusions of Theorem^ hold 

2.4.2 Thermalization and decoherence 

We describe some further results that could be of interest. The following result states that some observables tend to 
an asymptotic value as t — > 00. Let pg E be the (finite volume) Gibbs state at inverse temperature /3 corresponding 
to the interacting Hamiltonian H , i.e. 

1 

Let A £ £%(J%s) be a translation-invariant observable, that is, its kernel satisfies A(x' , e'; x, e) = A(x' + y, e'; x + y,e) 
for any y £ Z d (recall the notation introduced in Section [2.4. 1[) . For simplicity, we assume A to be an infinite- volume 
observable and whenever necessary, we make it into a finite-volume (A = Z d /L'Z d ) observable by restriction, i.e. 
A A := 1a^41a with 1a = ly <g> 1;2(A). Define the equilibrium expectation value of A as 

(A) = lim Tr[ Ps ^ A ® 1 E )] (2.32) 

Note that pg E depends on A, too. The existence of the limit follows easily by the expansions in Section [10. 31 

Theorem 3. Let A £ 38(J#s) be a translation-invariant particle observable, as described above. Under the conditions 
of Theorem^ (including the restriction on \X\), we have 

lim Tr[p g tA] = {A)b, where (3 := ft = (3 2 (2.33) 

t— >00 

If we drop the condition that fix = [3%, then the asymptotic state of the particle is given by a NESS (non- 
equilibrium steady state): 

7 



; e-f* H (2.31) 



Theorem 4. Let A be a translation-invariant particle observable. Under the conditions of Theorem^ we have 

lim Tiips.tA] = (A) ncss , (2.34) 

t— >oo 

for some linear functional A — > (A) ness on translation invariant observables, satisfying |(A) noss | < \\A\\ag^ s ^, 
(l) ncss - 1 and (A) ncss >0forA>0. 

The easiest way to give more details on the NESS consists in using the fact that it is a small perturbation (in 
A) of the NESS that one obtains in the Markovian approximation to our model, and this will be described below. 
Let us however quote one important property of both functionals Oncss and (-)p, namely decoherence. 

Choose the observables A y with kernel A y (x' , e' , x, e) = 5 x i- X)V , then there is a 70 > 

(A y ) ncss < Ce-nM (2.35) 

for some constant C, independent from y. The same statement holds for (A y )p as well, but this does not require 
our analysis since it is a property of the interacting Gibbs state introduced above. 

2.4.3 The Markov approximation 

As mentioned already above, we can describe our results qualitatively by referring to the Markov approximation 
to our model. Strictly speaking this Markov approximation retains the quantum nature of the problem, but as far 
as the the long-time properties of the system are concerned, we can describe what is happening with the help of 
the 'classical' Markov process that was already introduced before Assumption [Cl by specifying the rates j(-, •) on 
the state space T = <j(H sp i n ) x T d . One should think of the elements in F as 'good quantum numbers' for the 
Hamiltonian H§, with k £ T d momentum and e £ cr(iJ sp i n ) spin. Since we assumed it to be irreducible, this Markov 
process has a unique invariant state given by an absolutely continuous positive measure. We call the associated 
density HQ{k, e). 

If j3\ = $2, then we simply have 



/'Q 



(k, e) = (l/27r) d e-^ e (^ e"^')" 1 (2-36) 



The fact that this Gibbs state is uniform in k is due to the fact that the kinetic energy was chosen so small that 
on the time scale A -2 for which the Markov approximation is valid, the kinetic degrees of freedom do not couple 
directly to the phonons. 

For /?i t^/32, we do not have an explicit expression for the invariant density /xq(/c, e), but it is connected to the 
NESS discussed in the previous sections as follows. For a translation invariant A, we define 

(A) ncss , Q := 7i/ dfc E MQ(M]Te ifea; A(0,e, x,e) (2.37) 

T d e X 

Then we have 

(A) neas - (A) nesStQ = o(|A|°), A^O (2.38) 

In fact, the same statement holds true when /3± = /?2, but in that case it follows simply by explicitly comparing the (S- 
part of) the full interacting Gibbs state with fJ,Q(k, e), i.e. no analysis of the dynamics is necessary. Furthermore, we 
note that by the invariance property j(e,' dfc'; e, k) = j(e,' d(k'—k); e, 0), we can deduce that HQ(k, e) = (1/2tt) (j,Q(e) 
(independent of k). 

From the Markov approximation, we can also infer a diffusion constant. Assume that the Markov process on 
JF describes a particle that jumps in its (fc, e) coordinates and, in between jumps, propagates freely in space with 
velocity 

v % = 2m p 1 sinfc,, k £ T d (2.39) 



(note that this is the derivative of the dispersion relation m p l ^ i=1 (2 — 2cosfci) of A 2 -ffkin)- The diffusion constant 
Dq of such a particle is given by the velocity- velocity correlation function 



DgSij = ^ / dtiviWvjiO)) (2.40) 



where the expectation (•) is computed with respect to the stationary Markov process. Then, the diffusion constant 
D* in Theorems [U [2] has the asymptotics D* = X 2 Dq + o(A 2 ) as A — > 0. 

2.5 Related work 

2.5.1 Classical mechanics 

Diffusion has been established for the two-dimensional finite horizon billiard in [BJ. In that setup, a point particle 
travels in a periodic, planar array of fixed hard-core scatterers. The finite-horizon condition refers to the fact that 
the particle cannot move further than a fixed distance without hitting an obstacle. 

In [17], the hard-core scatterers are replaced by a planar lattice of attractive Coulombic potentials, i.e., the 
potential is V(x) = — X^ez 2 \ x - l • ^ n ^ a ^ case > the motion of the particle can be mapped to the free motion on a 
manifold with strictly negative curvature, and one can again prove diffusion. 

Recently, a different approach was taken in [3 : Interpreted freely, the model in [3] consists of a d = 3 lattice 
of confined particles that interact locally with chaotic maps such that the energy of the particles is preserved but 
their momenta are randomized. Neighboring particles can exchange energy via collisions and one proves diffusive 
behavior of the energy profile. 

2.5.2 Quantum mechanics for extended systems 

The earliest result for extended quantum systems that we are aware of, |20) . treats a quantum particle interacting 
with a time-dependent random potential that has no memory (the time-correlation function is 8(t)). Recently, this 
was generalized in |16j to the case of time-dependent random potentials where the time-dependence is given by a 
Markov process with a gap (hence, the free time-correlation function of the environment is exponentially decaying). 
In [21], a quantum particle interacting with independent heat reservoirs at each lattice site was treated. This model 
also has an exponentially decaying free reservoir time-correlation function and as such, it is very similar to [T6] . 
Notice also that, in spirit, the model with independent heat baths is comparable to the model of j3], but, in practice, 
it is easier since quantum mechanics is linear. 

The most serious shortcoming of these results (except for [8] , already discussed in the introduction) is the fact 
that the assumption of exponential decay of the correlation function in time is unrealistic. In the model of the 
present paper, the space-time correlation function, £(x, t), is the correlation function of freely- evolving excitations 
in the reservoir, created by interaction with the particle. Since momentum is conserved locally, these excitations 
cannot decay exponentially in time t, uniformly in x. 

In the Anderson model, the analogue of the correlation function does not decay at all, since the potentials 
are fixed in time. Indeed, the Anderson model is different from our particle-environment model: diffusion is only 
expected to occur for small values of the coupling strength, whereas the particle gets trapped (Anderson localization) 
at large coupling. 

Finally, we mention a recent and exciting development: in [llj . the existence of a delocalized phase in three 
dimensions is proven for a supersymmetric model which is interpreted as a toy version of the Anderson model. 

2.5.3 Quantum mechanics for confined systems 

The theory of confined quantum systems, i.e., multi-level atoms, in contact with quasi-free thermal reservoirs has 
been intensively studied in the last decade, e.g. by [2J [TU [TO] [22J . In this setup, one proves approach to equilibrium 
for the multi- level atom. Although at first sight, this problem is different from ours (there is no analogue of diffusion) , 
the techniques are quite similar and we were mainly inspired by these results. However, an important difference is 
that, due to its confinement, the multi-level atom experiences a free reservoir correlation function with better decay 
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properties than that of our model (if one assume enough infrared regularity) and the 'Markov approximation' to 
the model has a gap, whereas in the extended system it is diffusive. 

2.5.4 Scaling limits 

Up to now, most of the rigorous results on diffusion starting from deterministic dynamics are formulated in a scaling 
limit. This means that one does not fix one dynamical system and study its behavior in the long-time limit, but, 
rather, one compares a family of dynamical systems at different times, as a certain parameter goes to 0. The precise 
definition of the scaling limit differs from model to model, but, in general, one scales time, space and the coupling 
strength (and possibly also the initial state) such that the Markovian approximation to the dynamics becomes exact. 
In the model of the present paper, one can do this by considering the dynamics for times t = 0(A~ 2 ) and then taking 
A — > 0. The result is a Markovian approximation which was already referred to in Section l2.4.3l Of course, the point 
of the present paper is that we go beyond the Markovian approximation and we describe the dynamics for infinite 
times at fixed A. There are quite some results on scaling limits in the literature, for example [24l [T4l fT3l IT9l Il8l IT2] . 
and we do not attempt an overview. 

2.6 Strategy of proof 

We will now give a road map for the proof of Theorems Q] and [2] It is based on a careful analysis of the long time 
properties of the reduced particle dynamics given by the operator Z t defined in eq. (|2.14p and extended to infinite 
volume in Section [5) We view Z t as a bounded map Z t : £$i(J%s) —>■ £%i(J#s) i.e. 

The main ingredient will be a proof that in a suitable norm the rescaled large time limit 

lim Sr t Z t (2.41) 

existqj. Here S n is an operator implementing the scaling of the particle position occurring in Theorem [T] in the 
space 38{3§{3#s)) explained in detail in Section |3"U1 This large time limit is controlled in two steps. 

2.6.1 Random walk in random environment 

The first step consists of studying the dynamics on the time scale 0(1/ A 2 ). This scale is large enough so that the 
dissipative effects that we want to exhibit are clearly visible, and small enough such that a simple Duhamel expansion 
can be controlled. Thus, let us fix to = A _2 to with to of 0(1). The evolution of the system plus environment up to 
this time is given by 

U:=c-' ltoL . (2.42) 

Denote the reduced time evolution on this scale by 

T:=Z t0 . (2.43) 

In Section fl 1.1 1 we show that for A small enough, 

T = e t (-iL s +A 2 A/) + o{\\\ 2a ) (2.44) 

in an appropriate norm, cfr. Proposition 1 1 1 . ll Here, Lg = ad(iJg) and M is the generator of a "quantum Markov 
process" that is very closely related to the Markov process discussed in Section 12.41 

We will now compare the full evolution (|2.42[) to the one where the particle dynamics is given by the reduced 
evolution T and the environment evolves freely: 

U = T ® e - itoLE + B (2.45) 



3 This statement is of course only true in infinite volume 
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where L E = ad(_ff E )- This defines the "excitation operator B" acting on 3§i(J%s <8> <%?e)- This operator is analyzed 
in Section [TUJ There we show that B is small and "weakly correlated" (Propositions I10.T1 and fltT!?)) . To explain 
what we mean by this, consider the full evolution for times t longer than to- Taking t = Nto we have 

e -iNt L = U N = p g e -it L E + gjJV ^.46) 

We rewrite this 

u n = e -im L E( j + B(N))(T + B(N - 1)) . . . (T + B(l)) (2.47) 

where we use the shorthand T for T ® 1 and e _lt ° LE for 1 ® e ~ ItoLE and define 

B( T ) ■- e ^toL EBc -i(r-l)t L E ^ ^.48) 

Using cyclicity of the trace and the fact e~ lto E p^ = pg , i.e. the invariance of the environment state under the 
uncoupled dynamics we get 

Zm ps = Tr E [(T + B(JV)) . . . (T + B(2))(T + B(l))(p s ® p^)] (2.49) 

If B(t),t = 1, . . . , N were set to zero in (|2.49[) . the reduced evolution would be given by Z^to — T N i.e. by a 
discrete-time semigroup acting in the system space. This is our 'quantum random walk'. Using (|2.44[) . it is easy 
to show that it is diffusive; this is done in Section [TT] More precisely, we show that the process generated by M is 
diffusive and the diffusivity of T N follows then by simple perturbation theory. The presence of B{t) produces time 
dependence and dependence on the environment variables that are traced over in the end. We wish to think about 
the latter as time dependent noise and the Tr E as an expectation over the noise. Thus, given D £ 3S{38\(3#s®<ffiE)) 
we define E(D) g ^(^ X (J^ S )) by 

E(D) PS := Tr E [D(p s ® p r E of )] . (2.50) 

Using this notation we have 

together with 

and 

E B(t) = (2.53) 

for all r which follows again from the invariance of pg f . In Section |3. II we generalize this expectation to a larger 
algebra needed to analyze (J2.49I) . Given a set A = {t\,T2, ■ ■ ■ , T m } C {1,2,..., N} with the convention that n < 7"j+i 
we define the time-ordered correlation function 

G A := E (B( Tm ) ®s B( Tm _i) ®s • ■ • ®s B(n)) € [^(^(JT s ))]^ m (2.54) 

Here (gig denotes a product that is tensor product in the system space and operator product in the environment 
space, see Section 13.11 Like in classical probability we introduce in Section 13.21 connected correlation functions or 
cumulants G C A in our non-commutative setup. Propositions I10.T1 and fTO.21 show that these cumulants are small and 
inherit decay in time from the environment correlation function. 

Thus we want to think about T + B(t) as transition probability kernels in the system space which are random 
due to the environment dependence and we want to prove average or "annealed" diffusion for the associated process, 
a quantum version of random walk in random environment. 

2.6.2 Renormalization group 

The second step in the control of the large time asymptotics of Z t is to control the large N asymptotics of (|2.51|) 
composed with the scaling. This is achieved using Renormalization Group (RG) method which consists of studying 
(|2.49p in an inductive way. Pick an integer £ and define the RG map 

R{U):=S(,{U e ), (2.55) 
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Z NtQ =EU N (2.51) 



T = EU (2.52) 



where the scaling Se = Si (g)l acts on 3§(3§i(J#s)) only (see Section [373]) . Iterating it n times, and using (S() n = S^ 
we get 

U n :=R n (U) = S t ,{U l2n ). (2.56) 

Define now 

T„:=E[W n ]. (2.57) 

By (|2.5ip T„ gives the (rescaled) reduced dynamics at time £ 2n t : 

T n = SlnZf2n to . (2.58) 

Thus the existence of the limit f|2.41[) is closely related to proving (in a sense to be made precise) 

lim T n = T* (2.59) 

This is the content of Proposition 16.31 

T n is analyzed inductively in n. We define the renormalized noise B n in analogy with (J2.45I) by 

U n = T n ® e - if2 " toLE + B n . (2.60) 

B n is studied through its cumulants G c n A defined as for B. The RG map (|2.55[) leads to a recursion 

T n+l =F(T n ,G^.) (2.61) 

where the map F is given explicitly in eq. (|3.8p and another one for the correlation functions 

G n+i,A = ^A(T n , G c nm ) (2.62) 

with Fa given in eq. (|3.34l) . The recursion (I2.61|) is studied in Section [7] and (I2.62p in Sections [5] and [5] For this 
analysis a choice of suitable norm for G c n A is crucial: Norms which respect the structure of the iteration (|2.62p are 
constructed in Section [572] 

The convergence (|2.59[) is a consequence of the fact that that the map (|2.62[) contracts the norms of the correlation 
functions G c n # and preserve their temporal decay i.e. the noise is "irrelevant" . For the contraction one needs to 
study carefully the linearization of the map Fa (Section [5]). A crucial input to this analysis is the existence of two 
symmetries in the problem: unitarity and reversibility (Sections 14.11 and I4.3[) that provide the contractive factors. 

3 Renormalization group formalism 

In this section, we define the correlation functions G a and derive the recursion relations for their renormalized 
versions G n ^A as well as for the T n . We also discuss the symmetries preserved by the RG map. Throughout this 
section and the next one, all expressions are in finite volume. In particular, all sums that appear are finite. The 
infinite volume limit and the definitions of the tensor product spaces in that limit will be discussed in Section [5] 

3.1 Correlation functions of excitations 

We will now define the correlation function (|2.54[) and explain how the reduced dynamics (|2.49p can be expressed 
in terms of it. We abbreviate 

^s = #(#i (.*&)), ^E - #(#l(J&)) (3.1) 

such that U, B(t) are elements of Ms <S) £%e and T is an element of 38$. Define, for D, D' £ &$ <gi ^ E the object 

D <E) S D' £ 3? s <E) 3? s ® ^ E 

as an operator product in E-part and tensor product in S-part. Concretely, let first D = Ds<£)D-e and D' = D' S <E)D' E . 
Then 

D ® s D' := D S ®D' S ® D E D' E . 
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Extend then by linearity to the whole space S?s ® S&e ■ Iterating this construction we define for Di £ S?s ® S?e, 
i = 1, . . . ,m 

D m ®s • • • «>s D 2 ®s £>i e (^s)®" 1 <8> ^E- 
We define the 'expectation' 

as the partial trace 

E(F)(p s ,m <8> • • • <8 /Os,i) := Tr E (F(p s ,m <8> . • . <8> Ps,i <8> pif))- 

Hence we have explained the definition (|2.54[) of G^ . 

Note that by ([2351) G {r} = 0, and G A = G A +r since e - " ^^ = pg* . It will be convenient to label the S?ss 
and to drop the subscript S (since we will rarely need S?e) Let S$ T , r e N be copies of S?s- For A = {n, . . . , r m } 
with Tj < Tj+i, we define ,^4 by 

^U := ^ Tm ® f rm _, ® • • • ® ^ Tl . 

One sees that i?U is naturally isomorphic to S?® by identifying the right- most tensor factor to 3&\, the next one 
to ^2) etc. We use this correspondence to consistently view Ga as an element of Si a in what follows. Consider 
a collection A of disjoint subsets of N, then each of the spaces Si AeA 1S naturally embedded into SisuppA where 
Supp.4 = Ua^aA. Consequently, given a collection of operators Ka £ S?a, we can define 

<E> K A € Sg SuppA (3.2) 

AeA 

To express ()2.49p in terms of the correlation functions (|2.54[) we need one more operation, the contraction 
i : Si <E> Si — > Si. On elementary tensors define 

i(V®V'):=W (3.3) 

and then extend linearly. We can extend 6 to a map T : Si® — > Si in the obvious way: 

T(V m ® V m - X . . . Vx) := V m V m ^ ...Vi (3.4) 

i.e. T = L m ,m-i ■ ■ ■ *-3, 2^2,1 where Lk+i,k '■ Si® — > Si® equals i ® l^k-i . With these definitions eq. (|2.49[) can 
be written as 

Z Nt0 =T(E [(T + B(N)) ® S ---{T + 5(2)) ® s (T + B(l))]) (3.5) 

and expanding the product we end up with desired formula 

Z Nt0 = V T[ <8> T(t)®G a ], A c = {1,...,N}\A (3.6) 

* — ' tGA c 

AC{1,...,W} 

where the operators T(r) are copies of T £ Si embedded in Si T . Note that, for each operator appearing in the 
product, we have specified the space Si A or Si T in which it acts (indeed, recall that Ga was defined to act in Si A ). 
In contrast, the order in which we write the factors inside the T[-] does not have any meaning. For a completely 
explicit expression of T (and of Ta' , which will be introduced later) , we refer to Section 13.5.11 

3.2 Connected correlation functions 

The 'connected correlation functions" or 'cumulants', denoted by G C A , are defined to be operators on Si A satisfying 

G A ,= £ {aIa 01 ) (3J) 

partitions A of A' 

The tensor product in this formula is well-defined since Si A i = ® Si A whenever A is a partition of A' . Note that 

AeA 
this definition of connected correlation functions reduces to the usual probabilistic definition when all operators that 

13 



appear are numbers and the tensor product can be replaced by multiplication. Just as in the probabilistic case, the 
relations Q3.7P for all sets A! fix the operators G C A uniquely since the formula p.7[) can be inverted. 
Inserting (|3.7I) into (|3.6[) we get 



z Nt0 = T N + Vrff ® rw^fgG^) 



Ae<B{i) 



(3.8) 



where I = {1, 2, . . . , iV}, 05(7) is the set of non-empty collections A of disjoint subsets A of I, and Supp.4 = LUg^A 
Formula (I3.8[) follows from (|3.6[) by substituting ()3.7|> since, obviously, any disjoint collection A is a partition of 
SuppA The term T N in (|3TB|) originates from A = in (J3~B)) . 

3.3 Kernels and rescaling of space 

To define the scaling operator introduced in (J2.41I) we need to write our operators in a suitable basis. First, let us 
fix a basis for the space &i(J#s) . Since we assumed Jispin to be non-degenerate, we may label a basis of 5? by 
eigenvalues e € <r(H sp i n ). Hence ip € Mk = 5? ^ 2 (Al) may be identified with a function tp(x,e). Next, in this 
basis /o € 3§i{Ji(k) becomes a kernel (matrix) 

P = p(xL,e L ;x R ,e R ), 

where xl,xr € A^ and eL,e R G cr(-ff S pin) (R/L for "right/left"). The scaling operator will act on a particular 
combination of Xl and xr. Therefore we need to introduce new coordinates. Note that the construction below 
depends on the finite volume A (by the parameter L) in a trivial way, and we will not indicate this dependence 
explicitly. 

Let, for a vector a £ 2^ d > L a J s t an d for the vector with components [a\i = [a^J (largest integer smaller than 
Oj). Then we define 

,x L +x R ,xl-xr, 
x-=[ ^ J' v:= ^ 2 * ^ ' 

and 

7? := x L - (x + v) G {0, l} d (3.10) 

such that we have 

xl = x + v + rj xr — x — v (3-H) 

The variables eL,e^,ry will play a minor role in our analysis. Together with v, we gather them in one symbol s 
which takes values in the set 

S := a(H spin ) x a{H spin ) x {0, l} d x Z d /LZ d . 

The new variables (x, s) thus run through the set 

A := Z d /LZ d x S, (3.12) 

and we may identify 38(,3%s) with £°°(A ) i.e. p € ^i(J^g) becomes a function 

p = p(x,s). 
Similarly if € S% is given by a function on Ao x Ao i.e. a kernel K (x' , s'; x, s) which acts on 3§\{3%§) by 

(Kp)(x',s') = ^2 K(x',s',x,s)p(x,s) (3.13) 

(ap,s)eA 

We define a distinguished subset Sq C S 

Sq := {s = (e L ,e Rl rj, v) <E S \e L = e Rl v = 0, i] = 0} (3-14) 
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Note that 1S0 is independent of L, \Sq\ = dim ^, and that 

Trp = ^p(x,s)l 5o (s) (3.15) 

The scaling in the renormalization step affects only the variable x. At first, x takes values in Z rf /LZ d , but later 
on it lives on finer lattices. This motivates the definitions 

X„ := r n {Z d /LZ d ), K :=X„x5 (3.16) 

To reconstruct the original coordinates xl,xr from the (x, s) at scale n, one finds 

x L =£ n x + v + 7] x R = £ n x-v. (3.17) 

The scaling map transforms functions on A n into functions on A„+i, i.e. we define 

Si : l°°(A n ) -)■ ;°°(A„ +1 ) (3.18) 

with 

Sip(x,a) = £ d p(£x,s). 

Let us denote by M n the space of kernels K on A„ x A„ and the action 

Si : ffl n — > &n+l 

by SiK = S t KS~ l i.e. 

(S e K){x', s': x, s) = e d K(£x', s'; £x, s). (3.19) 

When summing over a variable x € X„, it is natural to take account of the cell volume i.e. consider Riemann sums. 
We do this by defining 

f dxf(x) := J2 t- nd m (3-20) 

Also, we use the shorthand z = (x, s) and the notation 



[ dzf(z)= [ dxdsf(x,s):= V r nd f(x,s). 



(i,s)eA„ 

These conventions are motivated by the fact that they preserve normalization. Thus we have 



(3.21) 



da; (S e f)(x) = / dx f(x). (3.22) 

and 

Trp = Tr5 f £p= / dxdsp(x, s)ls (s). (3.23) 

JA n 

Similarly the sum in (|3.13j) becomes an integral and the product of two kernels K\ , K 2 € £% n becomes 

(K 2 K 1 )(z',z)= [ dz"K 2 (z',z")K 1 (z",z) (3.24) 



We will also need a discrete delta-function 

5(z',z) = r nd 5 x , tX 5 s ,, a , z = (x,s),z' = (x',s') (3.25) 

which satisfies 

dz'5(z',z) = 1. 
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3.4 RG recursion for T 

Let us now study the renormalized reduced dynamics T n denned in (J2.57I) . T n belongs to the rescaled space 2% n . We 
copy the setup of Section [3TT1 defining spaces £% n ,A as products oi0l njT , copies of ffl n . We also define, analogously 

to (gag) , 

B n (r) :=c ^e 2rl toL EBne -i(r-i)f^t L E ^^ 

and the time-ordered correlation function 

G n , A :=E(5 n (r m ) ® s B n ( Tm -i) ®s • • • ®s B„(n)) € (#„)® m , (3.27) 

for A = {n, . . . , r m } with n < T2 < . . . < T m . With these notations, the concepts of Section l3~Tl correspond to 
n = 0, i.e. T and B(t) introduced there will now be referred to as To, Bq(t). As in Section [3. 1[ we embed G n ^ A into 
M n ,A and we write T n {r) to denote a copy of T n embedded in S% u ,t- 

The recursion relation for T n is obtained from the analogue of eq. (|2.51|) : 

T n+1 = S{T [E ((T„ + B n (£ 2 )) ®s ■ • • (T n + B B (2)) ® s (T n + B„(l)))] (3.28) 

which leads as in (|3.6p to 

T„+i= 51 S ^ r [ ® T «( T )® G n,A] (3.29) 

AC{1,...,£*} T€AC 

To get the analogue of p. 81 ) define first I T ', the set of times at scale n associated to a time r' at scale n+ 1, by 

7 T - := {£ 2 (r' - 1) + 1, ^ 2 (r' - 1) + 2, ... , *V} (3.30) 

Then 



T n+1 = S e [T r (]+ £ S ? T \( ® Tn(T))®(®Gl A ) 



Ae<s(i T , 



(3.31) 



where *B(I T i) is the set of non-empty collections A of disjoint subsets A of I T > , and Supp„4 = UaeaA. Note that 
the t' on the RHS is arbitrary. 

3.5 Recursion relations for correlation functions 

We derive the recursion for G^ +1 A ,, the cumulants at scale n + 1, in terms of those at scale n. The set of times at 
scale n that contribute to G^ l+1 A , is 

I A' :={J T , eA >lT> (3.32) 

We will need an extension of the contraction operator T. Given t' £ N, let 

IT' • *™Tl,I-.t ' ^n.T 1 



be the contraction as defined in Section I3T1 (followed by the imbedding of 8% n to M n ,T')- Given a set A' C N, we 
set 

Ta> ■= <S> T T > : M n j A , -> &„„# 

t'EA' 

Iii words, we contract within the time intervals I T ' ■ Below we illustrate the action of Ta> by defining it explicitly 
with kernels. 

Again, let *8(I A >) be the set of non-empty collections A of disjoint subsets A of Ia> ■ Then A £ Q3(/a') induces 
a graph Q A i (A) on the vertex set A' by the prescription that we connect t'^t'^ £ A' by an edge whenever there is 
an A E A such that 

An/ T ;^ 0, and A n J T / 7^ 0, (3.33) 

This leads to the basic relation between scale n and n + 1 



(» G^. A (» T„(r) 

Ae.4 r^Supp.4 



(3.34) 



G n+1,A' - 2^ <S^7a' 

.Ae*B( I A , ):Q A i (A) connected 

Here, the rescaling Si acts on all copies of & n , that is, we write Si instead of Si <g> . . . <E> Si. We omit the purely 
combinatorial proof of this relation. 
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T A 



<g> K A 

AeA 



3.5.1 Kernel representation and the map Ta' 

In Section [3. 3i we introduced coordinates x £ X„,s e5,z£ A„, such that operators in M n can be identified with 
kernels on A„ x A n . In a similar way, we can identify operators in !% n ,A with kernels on K n ^A x &n.A, where A KiJ 4 
is the cartesian product x T gjA„ iT and A„, r are copies of A„. Throughout the paper, we denote by za, z' a elements 
of A„, tA , and, hence, for any Ka £ &n,A, we have a kernel K A (z' A , z A ). 

Consider a partition A of A' (hence in particular A £ <8(/a')) an d consider a family of kernels K A (z' A , za) with 
AeA. Let us call a set of two consecutive times {t, t + 1} £ Ia' a bulk bond whenever {r, r + 1} c I T ' for some 
t'. Then we have 

(z' A „z AI ) = [dz' lA ,dz lA , [] S(z t+1 ,z' t )1[Ka(z'a,za) (3.35) 

{t,t + 1} bulk bonds Ae^l 

II <^T'> Z max/ T ,)<^T',ZminJ T ,) (3.36) 

where the discrete 6- function was defined in () 3 . 2 5 [) . 

4 Unitarity and Reversibility 

The correlation functions G c n A satisfy identities that are crucial for the proof of our results. The first of these "Ward 
identities" (called the 'unitarity' identity below) is due to invariance of the trace Tr(-) on 9S\{M'') under unitary 
maps and time-independence on the reference state of the environment w.r.t. the uncoupled dynamics. The second 
one (the 'Gibbsian' identity) is due to the invariance of the equilibrium Gibbs state under the dynamics. To discuss 
the latter we need to elaborate on the relationship between free and interacting Gibbs states, see Section H.2.11 The 
Ward identities will produce some additional smallness in our expansion. The 'Unitarity' identity eliminates certain 
terms exactly, the 'Gibbsian' identity does this only approximatively, as long as the temperature is not infinite. 

4.1 Ward identity from unitarity 

We start by examining the consequences of the conservation of probability, i.e. unitarity of e -1 *^. Let pse G &$\{3ri?) 
and abbreviate p — Sgn-ipsE- Then 

TrS„(r)p = 0. (4.1) 

Indeed, from (|H0|1 . ([235]) and (|3~25|1 we infer 

B n (r)p = S^iUnpU- 1 ) - (T n ® l B )p 

where U n € SS{J^°) is unitary. Thus, by using (|3.23|) . 

TrB„(r)p = Tr(U n pU- 1 ) - Tr((T„ ® l E )p) 

= Tr/9-Trp = (4.2) 

By picking appropriate pse this identity leads to the following one for kernels: 

dz' r l So (s' T )G c nA {z' A ,z A ) =0, t = max A, z' T = (x' T ,s' T ). (4.3) 

We rewrite this relation in the form in which it will be used. Let P £ 3§(l°°(S)) be a one-dimensional projector 
with kernel 

P(s',s) = u(s')l So (s), (4.4) 

where /i(-) is some function that satisfies (to ensure that P is indeed a projector) the normalization condition: 

E M s ') = L ( 4 - 5 ) 

s'eSo 
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In the rest of the paper, rank-one operators like P will be written as P = \p)(ls \- Then the Ward identity can also 
be expressed simply as 

dx' T {P®l...®l)Gl A (x' A ,x A ) = 0, T = maxyl (4.6) 



where P acts on (the s-coordinates of) 



^r-maxj 



4.2 Equilibrium setup: the case (3\ = /% 

In this section, we assume that the reference state of the two environments is in equilibrium at temperature /3 = 

01=02. 

4.2.1 Equilibrium states 

We already introduced the state /?£? on &§(Jffe). We now introduce a corresponding finite- volume reference state 
on £S(3%s). It is given by 

pf = (Tr[e-^])- 1 e -^s (4.7) 

and we define 

Ps 4 = pf^p^ (4.8) 

We will also need the Gibbs state of the coupled system. We introduce 

/4 = (Tr[e-nrV^, ^:=Tr E/ 4 (4.9) 

Hence we adopt the convention of using the superscript "ref" for states that are Gibbs with respect to the uncoupled 
Hamiltonian, and the superscript (3 for interacting Gibbs states and restrictions of those. In particular, note that 
Pg is not a Gibbs state but merely a restriction of a Gibbs state. 

To ease the discussion of the thermodynamic limit, we also introduce 

z/ cf := \A\pf, v* := |A|pg (4.10) 

As a consequence of translation invariance, the kernels of v° , z/' are constant in the variable x G Xo, hence they 
reduce to functions of s £ S; 

u p (x,s)=n l3 {s), if ei (x,s)=p, le{ {s) (4.11) 

Moreover, the chosen normalization ensures that 

J2 Ma) = J2 ^ rcf ( s ) = 1 ( 4 - 12 ) 

s£5o s£<So 

Next, we introduce the Radon-Nikodym derivative as an unbounded operator on Jff 

D rd := e |AF( / 3) e -4ff e |(if s+ H E ) ) (4 13) 

where Ai r (/3) is a number defined by 



e 



AF(/3) 



Tr [e 



-P(H S +H E )] 



Tr[e-^] (414) 

One can check that the 'free energy difference' AF((3) has a limit as A /• Z d (it is not proportional to |A|, because 
the number of particles does not grow with the volume). D r( j is an unbounded operator, even in finite volume, but 
it is obviously well-defined by the functional calculus. The reason why we call D r( j a 'Radon-Nikodym derivative' is 
that, formally, 

Pse = ^d/s^d (4-15) 



4.2.2 Correlation functions involving r = 

Let us define the operator Uq{0) by 

W (0)V := D Td ipD* d (4.16) 

such that formally 

pl E =Uo(0)(?£l (4.17) 

It is not hard to see that Uq(0) is densely defined on 8§i(,^f). This is however not our point. Rather we want to 
treat the operator Uq{Q) en par with the operator Uq introduced previously, hence we write 

W (0) = T o (0) ® 1 E + B (0) (4.18) 

where 

T (0):=EWo(0) (4.19) 

We can now extend formally the definition (|2.54p such that the set A can include t = 0. We define an additional 
copy, 3% n fi, of S£ n such that again G„ ^ is embedded into & n ,A- When defining the renormalization transformation 
of Wo(0), we omit the iteration and keep the rescaling only, such that we have 

U n {{)) := StnUo{0), B n (0) := S*.B o (0), T n (0) := S*»r (0) (4.20) 

With these definitions, we can also formally define the correlation functions G n .A € &n,A at higher scales n > 1. 

Lemma 4.1. TTie operators T„(0) and Gn^jG^^ wzi/i A 9 0, are bounded. The relation (J3.34I) remains valid for 
A3 if one defines Iq := {0}. 

The boundedness will be easily proven with help of expansions in Section 110.31 The validity of the recursion 
relation is formally obvious. For later use, we also define the bounded map Z t : 8$i(,3tf$) — > 3$i(J%s) 

Z tPs := TY E [c- [tL (D ldPs p^D; d )} = Ti- E [e- itL U a (0)( Ps ® pjf)] (4.21) 

4.3 Ward identity from equilibrium 

To state this Ward identity, it is convenient to introduce an additional piece of notation. Given K € 8% ® 8% and 
i/i G ^i(^l), let Kip G ^® @i(Jfs) be given by 

(V 2 <8> W := F 2 (8 (V x ^) (4.22) 

on elementary tensors K — V^ <S> V\ and then extending by linearity. This notation is only used in the lemma below 
and its proof. 

Let us also define the scaled versions of the states (|4.10p 

vf := S e S ct , i£:=W (4.23) 

Then we have the Ward identity: 

Lemma 4.2. 

E(B B (r) ®s B„(1))k£ = -E(B n (T)®s-B n (l)S n (0)X f -E(S n (T)®s(T„®l E )B n (0)K ef (4.24) 

+ E(B rl (r-l)® s B n (0)X of (4-25) 

Proof. Up to trivial rescaling, the proof is the same for all n, we will therefore for simplicity set n — and omit it 
everywhere. We also abuse notation by writing T„,T n (0) for T n <g> 1e,T„(0) ® 1e- We introduce the projector P 
that acts on 8§x(,^if) by 

■PV = (Tr E ip) (g> p™ { , ^ € #i(Jf) (4.26) 
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We start from 

B{l)V P ^ = -B(l)(l-7>)/4- 7^4 + 0^/4 (4.27) 

where we used B(l) = e ltoLE U — T and pg E = Wpg E . Next, we substitute 

PL = ^(O)Pse = W(0)(p r s of ® p r E cf ) = (T(0) + 5(0)) (pf ® p r E cf ) 

and use Ppgg = Pg <£> p E f to obtain 

S(l)(p^p r E ef ) = (- J B(l)(l-P) J B(0)-T J B(0) + c i *° LE J B(0))(p s cf ®pr) 

+ (e i ' oLE T(0)-rT(0))(p r s of «)p r E cf ) (4.28) 

This relation has the form Y^ Dj(ips,j <8> p E f ) = 0, for some ■i/>sj G 3§i(J#s) and £)j operators on ^i(J^), hence it 
implies that 

J2 E ( D ' ®S £>i)^S j = 0, for any D' (4.29) 

3 

where we used the convention (|4.22p . We choose D' = B(t) and we spell out (|4.29|) . obtaining 
E(B(r)(8)sB(l))pf = -E(B(t) ®s B(l)B(0))pf 

-E(B(t) <g> s TB(0))p r s cf + E(B(t - 1) <g> s B(0))pf (4.30) 

where the two first terms in (|4.28[) containing no -B(-)-operator have disappeared because E(B({t})) = 0, and we 
used E(B(r)e ltoLE (gig . . .) = E(B(r — 1) ®g . . .). The lemma follows upon multiplying with |A| and reinstating the 
subscript n. 

U 

4.3.1 Ward identity in terms of correlation functions G An 

As it stands, Lemma 14.21 is not written in terms of the correlation functions G C A . We can however easily rewrite 
it in that way. Since the Ward identity compares operators in different copies of & n , it is unnatural to index 
the z, z'-coordinates by r £ A and hence we use arbitrary indices, with the convention that the z, z'-coordinates 
corresponding to earlier times stand to the left of the later ones. We again drop the index n. 
The Ward identity reads 



&z a G c {lT} {z' a ,z b ';z a ,z b )vP{z a ) = j dz a L r (z' a , z b ' '; z a , Zb)v m (z a ) (4-31) 

with 

^(^1 z b ; z a , Zb) '■= — I dz G| 0;1]T j(z, z al z b ; z a , z, z b ) 

& z T{z a , z)G{ r j(z, z 6 ; z a , z b ) + G^ T _ 1 y(z al z b ; z a , z b ) (4.32) 



The verification of this relation from Lemma T4. 2 1 is by inspection, using the fact that Ga = G C A whenever \A\ < 3, 
which in turn follows from the vanishing of Ga whenever \A\ = 1. 

Remark. The reader should think of eq. (J4.31I) as analogous to the unitarity Ward identity, eq. (I4.3[) . Namely, the 
RHS turns out to run down to zero fast under the RG (i.e. asn-> oo) and this lets us get extra contraction for the 
RG flow of the two-point function. 
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5 The infinite volume setup 



Up to this point, our whole treatment was restricted to finite volume TL d jL7L d . This allowed us to neglect all sorts of 
analytical questions since all the operators pertaining to the S-part were actually finite matrices (this is not true for 
operators on J^e, though). In this section, we indicate which quantities remain meaningful in the thermodynamic 
limit and we give their precise definition. Note that the lattices X n ,A„ remain meaningful with L = oo, as was 
already indicated in Section l3~3l However, the spaces & n ,A, defined as tensor products in finite volume, are no 
longer uniquely defined and we will choose a norm to specify them in Section 15.21 

5.1 Thermodynamic limits 

Let us now indicate explicitly the volume dependence in the evolution operators T„ , T„ (0), the correlation functions 
Gn A^n A an d the 'unnormalized states' i^' A , fjj ' . Those objects were previously simply denoted by the same 
symbols without A but now we reserve this notation for infinite-volume quantities. The thermodynamic limit is 
then defined by point- wise convergence of kernels (well-defined since A^ is naturally a subset of A„) Recall the 
assumption ()2.17|) . it leads to 

Lemma 5.1. The kernels of the operators T„,T„(0), the correlation functions G^ A ,G^ A and the 'unnormalized 
states' i/£' , i/j™' converge pointwise os A/ 7L d . The limiting kernels satisfy the recursion relations (|3.31[) . (|3.34|) 
and the Ward identities (|4.3j) and (|4.31j) . In particular, the sums on A implicit in the RHS of these recursion 
relations are absolutely convergent. 

This lemma follows directly from the bounds of Lemma 110.11 and (for r = 0) those of Section 110.31 

5.2 Norm on kernels 

5.2.1 Tensor product completions 

We have introduced the space 3% n (on scale n) for operators acting on (rescaled) particle density matrices. As long 
as the volume A is finite, this space is finite-dimensional and hence isomorphic to &(3$i(J%s)). Similarly, the tensor 
product spaces 8%A,n are finite dimensional and thus all norms are equivalent. However, as A /* X d , care is needed. 
Even if we choose a norm for & n there are several ways to complete the algebraic tensor products. Our choice is 
dictated by the recursion relation (|3.34[) which we want to be continuous in the norm. The norms are of course 
meaningful for A finite as well, and in the definition below we will not distinguish different values of A. 

Recall that an element of 8% n ,A is represented by the kernel Ka = Ka{z' a , Za) with A — {t\, . . . , r m }, t% < tj+i 
and za = (2-7-1 > ■ • ■ 7 z r m ) G &n,A- Recall also that z T = (x T , s T ). Since the scale subscript n does not play any role as 
far as the definition of the norm is concerned, we suppress it in what follows and we will simply write ^^a, etc... 

We need to distinguish in the norm the x and the s variables and so will view ffi as the tensor product M x ® Stf/ 
where K 6 S% x is represented by the kernel K(x',x) and K £ M s by K(s',s). Introduce the following L 1 — L°° 
norm in the spaces M x and £# s 

11*11 :=m^(snpJdy(\K(y',y)l sup J &y'\K{y' ,y)\)\ (5.1) 

where y stands either for x if K £ S% x and for s if K G S% s . The space Ma can be viewed as the tensor product 
S%% m ®S% s Tm <£) ■ ■ .<g)& x 1 (S)&T 1 - It i s however simpler to consider more generally kernels K e (&?L 1 3%'' i where Ci € {x, s}, 
hence not necessarily with the same number of x- and s-coordinates (Note that we drop the r-labels since what 
follows is independent of the labelling of tensors). We now define the norm ||iif|| on such kernels, inductively in the 
degree m. 

Write K = K(y' ,y) with y — (yi, . . . , y m ) and yt equals Xi or Si depending on er^, and analogously for y'. For 
1 < i < m, fix y< = (j/i, . . . , 2/i-i) and y> = (y,+i, . . . , y m ) and let K l e & ai be the restriction of K to S% ai with 
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y< and y> kept fixed. Then Ki :— \\K l \\ is a kernel of degree m — 1. We define the norm inductively in the degree 
of kernels by setting 

|| jq := max piU. 

i 

The virtue of this norm is that it behaves well under tensorproducts and contractions, the basic building blocks 
of the the recursion relation Q3.34|) . The contraction analogous to (|3.3[) i : M^ <X> !%° — > M° is given by 

(t,K)(y',y) = fdyK(y,y';y,y) (5.2) 



We extend this map to ^L. x M ai in the natural way: for i,j s.t. <Ji = o~j let ty act in the «:th and the j:th factors 
M°* <g> ^" Tl . Finally, we write \K\ for the operator with kernel \K\(y',y) := \K(y',y)\. Points 1) and 2) of the 
following lemma are used throughout, point 3) is necessary only in the proof of Lemma 17.21 

Lemma 5.2. Let K and L be kernels as above. Then 

1) 

||ff®.£||<||.K'||||.£|| (5.3) 

2) Let oi = <jj such that tjj is defined, then 

\\ kj K\\ < \\K\\ (5.4) 

3) Let K be of degree m with o\ = o~2 = ■ ■ ■ = o~ m , then, for any i €E {1, . . . , m}, 

SUP |tl 2 l23 • • ■ i,n-l, m K\ < || SUp |/f||| (5.5) 

y>y' yi,y'i 

where on the RHS the norm is applied to a kernel of degree m — 1 . 

Proof. For the first claim, we proceed by induction in the degree of K ® L. Obviously (K ® L)i = JP <£> \L\ if 
i = degree(L) + j, and (K ® L) 4 = \K\ ® Li if i < degree(L). Thus 

||if<8.L||<max(||if i (g»L||, pr®J^||) < max(||^||||i||, ||^||||i,-||) = ||if ||||i|| 

where we used induction in the second step. 

For the second claim, the inductive definition of our norm means it suffices to check the claim for K of degree 
two. Thus consider the first term in (15. fl) for lK: 



sup / dy 

y' 



dyK(y,y';y,y)) 



< sup dy dy\K(y,y';y,y)\ < sup / dysup / dy\K(y",y';y,y)\ 

y' J J y' J y" J 



and the expression on the right is bounded by ||-Ki||. By analogous reasoning, the second term in (|5.1|) is bounded 
by H-ft^H and hence we have indeed \\cK\\ < \\K\\. 

To get the third claim, it is sufficient to check it for m = 2 and for m = 3 with i = 2, since all other cases can 
be reduced to one of these. We check explicitly the case m = 3 with i = 2. The LHS is estimated as 

sup dy dy\K(y,y,y';y,y,y))\ < sup / dy'" sup / dy sup \K(y"', y" , y'; y, y, y))\ 
y,y' J J y,v' J y J y" 

< sup / dy'" sup / dy sup \K(y"', y", y'; y, y, y))| 
v J y' J y,y" 

and one now verifies easily that the result is indeed smaller than the RHS. The m = 2 case follows by similar, though 
simpler reasoning. □ 



4 The contraction defined in (|3.3j) acts in both x and s coordinates, whereas the one defined in (|5.2j) acts in one 
of them 
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5.2.2 Definition of the norm || • || 7 

We will modify the norms introduced above to encode information about the decay properties of the kernels. Thus 
we will define a norm Hif^H^ that depends on two parameters: 70 (giving the decay in the v — v' coordinate, and 7 
given the decay in the x — x' coordinate. Note that the parameter 70 is not included in the notation for the norm, 
since we will never need to consider a different one. The parameter 7 will equal a multiple (of order 1) of 70. Given 
a kernel K A G &Ai let 

Kl(z' Al z A ) ■=K A {z A ,z A )e-t^eA\<-^\ e ioT. T eA\<-^\ ( 56 ) 

where we wrote z T = (x T , s T ) and s T — (v T , rj T , (ei,) T , (e/j) T ). Then we set 

\\K A \U:=\\Kl\\. (5.7) 

Note that 

(KaQKbY 1 =Kl®Kl. (5.8) 

Using triangle inequality \y' — y\ < \y' — y\ + \y — y\ in (15.21) we get for point-wise nonnegative kernels K A 

{Lrr>K A )~< < L TT .K\. (5.9) 

with the contraction i acting in both x and s coordinates (as defined in p.3p ). These observations lead to 

Lemma 5.3. 1) Let A be a collection of disjoint sets A C No such that I A > — Supp.4 (see Section \3.5\) and let 
for all A G A an operator K A £ S£ A be given, then 

(5.10) 
L^ J 7 IeA 

2) Let K A £,% A . Then 

(VI 1, 

and 

for 7 < i '. 

Proof. For 1) we bound (|5.10[) from above by replacing K A by \K A \ and then use (|5.8p and (|5.9p and the fact that 
T is a product of contractions ii 3 • . For the first claim in 2) observe that the norm (|5.1[) is invariant under Si and 
the effect in 7 is the stated one. The second claim is obvious. □ 

Remark From now on, the space £% is meant to be equipped with the norm || • || 7 . The exponential factors in our 
norm (provided 7 > 0) imply in particular that kernels K g Stf, map the space 

& p (Jf s ) = iv G #(J*&)| Tr(VV*) p/2 < 00} 

into itself, for any 1 < p < 00. 

5.2.3 Definition of the norm || • ||<y 

When studying phenomena like diffusion that are intimately connected with the translation invariance of the system, 
we will mainly need to consider the coordinates x,x' (see Section 15.3.11 on translation invariance). Therefore, we 
often want to consider K £ S% — S% x ® M s as a kernel in x, x' but taking values in the 'internal space' M s which for 
short will be called Sf . & is completed with the norm 



Ta> 


® K A 

_A£A 


< n 11^1 

T AGA 


\\S t [K A ]h = WKaW./z, 




\\KaW~, < \\Ka\W 



\F\\<,:= 



max ( supj] |F(s',s)|e 7o|,J - , '' l ,sup^ |F(s', s)\)e^ v - v '^ ) (5.12) 



23 



Obviously, the norm || ■ \\y relates well to || • || 7 on S%. Given K G 3?,, 

max ('sup I dx\\K(x\x)\\c$^ x - x '\sMV [ dx'\\K{x' ,x)\y G ^ x - x '\\ < \\K\\ 7 . (5.13) 

Indeed, looking back at Section IS .2. II and taking o\ = x and a-i — s, the LHS equals U-K^ll whereas the RHS equals 
max i= i !2 11-^711- 

5.3 Symmetries 

5.3.1 Translation invariance 

We implement spatial translations on Ms = 5? ® l 2 (Z d ) by the operators V u G ^(Ms), u G 1 d \ 

(V u iP)(y):=ilj(y + u)€y, i> e M s (5.14) 

Write in general Ad(S)0 = SOS^ 1 for 0,S E SS[M£). Then Ad(V^) implements translations on density matrices. 
In our Ao-coordinates, it acts as 

Ad(V u )p(x,s) =p(x + u, s) (5.15) 

We will also denote by Ad(Ki) the operator on A n defined by (|5. 15|) . now with x, u G X„. It is then natural to call 
a K G Stf, = S% n translation invariant if it commutes with translations, i.e. 

V U KV- U =K, for a e X n (5.16) 

For such a translation invariant operator K, we will often abbreviate the reduced kernel 

K(x'-x) = K(x',x) {K(x',x)€&) (5.17) 

and we have 

Lemma 5.4. Let K G 3% n be translation-invariant, then 

\\K\\ 1 = [ 6x\\K(x)ye lM (5.18) 



Proof. The inequality > follows from (|5. 13|) . Following the discussion after (|5.13[) . in order to get <, we need to 
bound ||i^7ll ^ 11-^2 II- This follows since the supremum over the x or x' coordinate can be dropped because of 
translation invariance and hence we get an upper bound by moving the supremum over the s or s' coordinate to 
the right. □ 

5.3.2 Fourier Transform 

The translation invariance suggests to Fourier transform the kernel in the variable x G X„. The dual space to the 
lattice X„ is the torus 

T„ := (£ n T) d (5.19) 

We define for p G T„ , 

p(p,s):= [ dxe ipx p(x,s), forp€^(A n ), (5.20) 



and, for K G & n , 

K(p):= I dxc ipx K(x), (K(p)€&). (5.21) 



where the sum on the RHS is absolutely convergent if || K\\ 7 is finite for some 7 > 0. It follows that, if we have two 
translation invariant kernels K\,Ki taking values in Sf, then 

Kjr 2 (p) = £i(p)# 2 (p) (5.22) 
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where the product on the LHS is in M n and on the RHS in Sf . The following standard consequence of exponential 
decay will be used throughout. If if is a translation-invariant kernel then 



sup \\K( V )\\y<\\K\\^ (5.23) 

Imp<7 



5.3.3 Symmetries of the lattice 



Let us investigate the action of a lattice symmetry O on our operators. We consider a density matrix p, and we 
abbreviate the transformed density matrix as po = Ad(Vb)p, such that 

Po{xl,x r ) = p{Ox Ll Ox R ) (5-24) 

To write this transformation in our new coordinates x, s, we let xo,vo,tio be the coordinates corresponding to 
Oxl,Oxr (we suppress the coordinates eL,en since they are untouched by O). Define the vector e by 



(5.25) 



a _ f if (0 V ), e {0, 1} 
*' ~ I 1 if (Ov)j = -1 
Then, one checks that 

xo = Ox — e, vo — Ov — e, Vo = Or] + 2e (5.26) 

Hence po( x , v , v) — p{Ox — e, Ov — e, Or] + 2e). We compute the fourier transform of p in the variable x. 

poip, v, V) = e ipe >(0-V, Ov -e,0 V + 2e) = (i^oPiO^p, ; ■)) («, »?) (5.27) 

where J P! o is an invertible transformation acting on the degrees of freedom v, r\. 



6 Convergence to a fixed point 



We now state the induction hypotheses for the RG flow of the reduced dynamics T„ (Proposition 16. 1[) . and the 
correlation functions G c n A (Proposition 16 . 2| ) . Their validity for all n implies our main results, as we show below in 
Section 16.21 The inductive proof of these induction hypotheses is postponed to Sections [Tl llll and we provide a brief 
guide to the proof in Section 16.1.31 

6.1 Induction hypotheses 

Let us first discuss constants and small parameters of our theory. The Hamiltonians Hs, He, Hi as well as the 
environment correlation function £ and temperatures /3, j3i are considered fixed and constants depending only on 
them are denoted by C (large constants) and c (small constants). The adjustable parameter in the problem is the 
coupling A and in the proof also the RG scaling factor £ and To entering the initial time scale to = A~ 2 to. 
We define the exponents 

J (2a - l)/4 1/2 < a < 1 J not defined 1/2 < a < 1 

"' [ (4a-l)/8 l/4<a<l/2 ' "* '" \ a/4 1/4 < a < 1/2 (< '""' 

with a the correlation decay exponent in Assumption El The exponent a\ (I for 'initial') appears in the treatment of 
correlation functions G C A with A 9 ('boundary' correlation functions; those were defined only in the case /?i = fc) 
and we define this exponent only for a < 1/2. Then, we introduce the 'running coupling constants' 

e n :=r na e , e :=C|A| Q (6.2) 

ei,„ = r nSl d l0 , ei, = C|A| 2 - 2q (6.3) 
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In the Sections [3 [8] and |9] the A appears only through 60,61,0 which can be considered the fundamental small 
parameters. 

Our basic convention is that £ is chosen large enough compared to the fixed parameters and we will for example 
freely assume that t~ 1 C < 1. The coupling constant |A| (or 60,61,0) is then chosen small compared to £, such that 
we can freely assume that |A|C(^) < 1 for quantities C(£) depending on £ but not on A. Finally, unless otherwise 
stated all constants are uniform in n, the RG iteration. 

6.1.1 Induction hypothesis for T n 

The first induction hypothesis concerns the reduced evolution T n and T n (r = 0) of (|4.20[) . The hypothesis involves 
parameters po, 7o and Do that will be fixed in Section ITT1 they result from the weak coupling analysis. 

Proposition 6.1 (Induction hypothesis for T n ). There exist tq, £q < 00 and X(£) > such that for £ > £q and 
|A| < X(£) the following holds, uniformly in to G [to, 2to] and in n. 

1) Analyticity . The operators T n (p) and T n (0,p) (X)' refers to r — 0) are analytic in a strip of width 70 with 

sup \\f n (p)\\ 9 <C, sup ||f n (0,p)||* < C (6.4) 

|Imp|<7o |Imp|<7o 

Moreover, T n (0,p) = T n -i(0,p/£) for n > 1. 

2) Diffusion. Let |Rep| < p„ given below and |Imp| < 70. The operators T n (p) have a simple eigenvalue 

e f ^ p \i.e. 

R n (p)f n (p)=e^R n (p) (6.5) 

where the one- dimensional spectral projector R n (p) is bounded as 

\\Rn(p)h < C. (6.6) 

and the complementary part as 

||(1 - Rn(p))f n (p)y < \£~> =: b n , (6.7) 

The analytic function f n (p) satisfies 

\f n (p)+D n p 2 \<ri n \p\ 3 (6.8) 

Moreover, for n > 1, 

\\R n {p) - R n ^(p/£)\\^ < V^rft"" 1 ), \ Dn _ Dnl \ < ^r £-i^~V (6.9) 

The constant p n is given by 

e -hD n pl =r %n e -±n oP l_ (610) 

and po < Dq/2. The spectral projection R n (0) is given by 

iZ„(0) = |/iT n ><lsol (6- 11 ) 

(notation as in Section ^. 1\ ) where ux n defines a probability measure on Sq: ^2 se g fJ-T n (s) = 1 and Ux n (s) > 0. 
Moreover in the equilibrium case fix = fi% = f3 

\\R n (0)-\^)(l So \\\ w <C^ (6.12) 

with pJ as in Section\4-2.1\ 



'For example, consider T n (p) as an operator on l°°(S), see the discussion in Appendix lAl 
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3) Gap. Let |Rep| > p n and |Imp| < 70. Then 

\\fn(p)\\ 9 < b n (6.13) 

4) In position space we have 

\\T n (x)\\ 9 < Ce-wl*! (6.14) 

For later use we note that ()6.8|) together with po < -Do/2 implies 

e -|D„(Rcp) 2 -C(Imp) 2 < | e /„(p)| < e -i£>„(R.cp) 2 +C(Imp) 2 _ ,g ^g-j 

6.1.2 Induction hypothesis for G^ A 

We move to the correlation functions G c n A . First, we define a distance-like function on sets A of times as follows: 
assume that A = {ti, . . . , r m } with ri < . . . < r m , then 

m— 1 

dist(A) = dist(n, . . . ,r m ) := JJ (1 + |r J+1 - tA) (6.16) 

and we will always use dist(A) Q = (dist(A)) Q to quantify the decay in time of the operators G c n A , with a as in 
Assumption |A] Since G c n A are translation invariant in time if ^ A it suffices to consider two cases: min A = 
and min A = 1 . 

Proposition 6.2 (Induction hypothesis for G C A ). Let to, £, X be as in Proposition \6. 1\ and recall e n defined in (16. 1|) . 

(a) Letl>a> 1/2. Then 

J2 dist(Ar||G^|| 107o <e^, fc>2 (6.17) 

AcN:| J 4|=fc,minA=l 

(b) If the environment is in equilibrium, i.e., if Pi = $%, then (a) holds for 1/4 < a < 1/2. In that case the 
correlation functions with Ab satisfy 

J2 dist(Ar\\G<; hA \\ Wl0 <e hn et k>2 (6.18) 

AcN:|j4|=fe,minA=0 

6.1.3 Plan of proof of induction hypotheses 

Given the induction hypotheses Proposition ^. II and Proposition 16 . 21 up to scale n, we prove Proposition 16. II on scale 
n + 1 in Section [7] This result is stated explicitly in Section 17.31 

The proof of Proposition 16. 21 on scale n + 1 is spread over Sections |8] and HI In the former section, we treat the 
linear part of the recursion relation and in the latter the nonlinear part (notions not yet defined) . The nonlinear 
part is straightforward and one does not need to distinguish between the cases a > 1/2 and a < 1/2. In particular, 
the equilibrium condition j3\ = p2 plays no role here. The linear part is slightly tricky; to treat the case a < 1/2, 
we have to exploit the equilibrium condition and consider correlation functions with A 3 0. The linear part also 
dictates the choice of the exponents a, a\ in eqs. (|6.2I6.3[) . respectively. 

Then, we need to establish the Induction hypotheses on scale n = 0. In Section [TU1 we establish Proposition 
16.21 and in Section QTJ we establish Proposition 16. II In those sections, we also outline how to choose the constants 
7o,p and D . 
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6.2 Proof of the main theorems 

We assume the Induction hypotheses Proposition 16. II and Proposition 16.21 for all n. Then we have 



Proposition 6.3. Let to,£, A be as in Proposition \6.1\ There is a projector P* e 'S of the form, \fj,*)(lg \, and a 
diffusion constant D* > such that for p in the strip \Imp\ < 70 

fnip) - e' 3 * p2 P* <C£- cn (6.19) 

for some exponent c > 0. In the case fi\ = /?2, we have p* = p" , the projected Gibbs state \4- 11\) . 

Proof. Let first |Rep| > p„. Then the claim follows from the definition (j6.10p and the bounds (|6.13|) and (|6.15[) . 
For |Rep| < p n note first that by (|6.9[) both R n (Q) and D n are convergent sequences. Calling their limits P*,D* 
respectively we get 

\\R n (0)-n\vAD n -D*\<el/ 4 . (6.20) 

Let m — en where < c < 1 is such that \£~ m ^ 2 p\ < 70. Such a c exists since p n > |Rep| and p ra grows not faster 
than (Cnlog^) 1 / 2 with n. Write 

m— 1 

\\Rn(0) ~ R n ( P )h < E \\Rn-U + l)(p£- U + 1) ) - Rn-oipi^h 

3=0 

+ ||i?„(0)-i?„_ m (0)||^ + ||i 1 , „_„ l (^-" l )-i?„_„ l (0)||^ (6.21) 

The first and second term on the RHS is bounded by e n _ m = e«_ % using (|6.9p . The third term is bounded by 
C(7o) -1 ^ m / 2 by analyticity of p — > R n - m (p) in the ball of radius 70 at origin. Hence 

\\Rn(0)~R n (p)\\w<Cr cn . (6.22) 

Next, bound 

< C(\f n (p) + D nP 2 \ + \D n -D*\)<Cr cn (6.24) 

by (|6.8j) and (|6.20|) . and the slow growth of p n . 

Combining (|6T20|l . (|6~22|) and (f6T24|) with (f6?f|) yields the bound (|6T9]) . It remains to argue that P* = P fj in the 
case /3i = /?2 = /3- Take A finite and recall that pg E = U n= o(r = 0)pgg and, by the invariance of the Gibbs state 

pi = Tr E [e-" L p^ E ] = Z t pf, t > (6.25) 

with Z t defined in Section 14.2.21 On the other hand, 

S/» [Z pnta ] = E(U n U n (0)) = T n T n (0) + T[G c m} J (6.26) 

where (0) refers to r = 0. On both sides, we can take the thermodynamic limit. Moreover, we know that 
||G? x i „||i07o < e_c ™ by Proposition 16.21 Let us then Fourier transform (16.26[) and apply the p — 0-component to 
p Icl , this yields 

/=f„(0)f„(0,0K cf + O(e- c ") (6.27) 

where the LHS follows from (jo^ and we wrote f„(0,0) = f n {p = 0,r = 0),f n (0) = f n (p = 0). In fact, by (j6~^5j) 
for t = 0, we have also f n (0,0)// of = // and by Proposition [Ol we have f„(0) = P* + 0(e _cn ). Hence 

/ = P*// (6.28) 

and this of course implies P* — P^ . 

U 
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6.2.1 Proof of results in Section 12.41 along a subsequence of times 



We argue that Proposition 16.31 implies the results of Section T2.4I along the sequence of times i„ :— £ 2n to. The 
resulting diffusion constant is given by 

D* = to l D* = t^X 2 !)* (6.29) 

Let us first write the claims of Theorems Q] and [2] in terms of the RG. Express the time i in units of the kinetic 
time scale, t = sio, to = A~ 2 to- Then we have 

Tr[e ip ^ s , t ] = Tr[e iP ^S^[Z sto ]^p s ,o] (6.30) 

For t = t n = t 2n to, S/jj[Z sto ] — T n and by the Fourier transform, we get 

Tr[e ip ^p s , tn ] = tr[f n (^=)p Sj0 (-^=)] (6.31) 

where the 'trace' tr is defined by tr[^] = ^2 s£S ip(s). Both T n (by Proposition 16.31) and ps.o (by the finite-range 
condition (12.111) ) are analytic in p and uniformly bounded in the strip |Imp| < 70. Hence, by Proposition 16. 31 (J6.31I) 
converges as n — > 00 to 

e- flV tr[P*Ai,o(0)] = e- DV . 
By the Vitali theorem, the derivatives converge as well and Theorem [T] and [2] follow (along a subsequence). 

For Theorems |3] and |4] we need to consider translation invariant observables A instead of e ip ^* in (J6.31I) . Its 
kernel A(x, s) = A(s) is constant in x and so we get from Proposition 16.31 

lim Tr Ap s e . n = ti[AP* p Sfi (0)} = tr[Ap*} (6.32) 

In the case where p* — pP , the last line of course equals (A) p. The decoherence result (|2.35p is a simple consequence 
of the fact that \\P*\\y < C, hence £ s p*{s)e^ v \ < C. 

6.2.2 Proof of results in Section 12.41 for general times 

Let us now consider general times in (J6.30I) . We will generalize a bit our RG iteration. It is easy to see that we can 
run the RG iteration i.e. Propositions 16.21 and 16.31 as well with scaling factor 2£ (by possibly reducing A) and we can 
also at each iteration step choose arbitrarily between the two factors. Let a € {0, 1} N label the possible choices: if 
a n = at the n:th step apply scaling £, if a n = 1 apply scaling It. Denote the resulting sequence of operators T nM . 
They satisfy Proposition 16.31 uniformly in a and hence we obtain projectors P* and diffusion constants D*. 



Lemma 6.1. P*,D* are independent of a. 

Proof. Denote the corresponding times by t n ^ a : 

t„,a = A- 2 £ 2 "4 m t (6.33) 

where m = m(n,o~) := #{i < n : o~i — 1}. Since T n „ = S h^t n we have Tn a — Tn a' whenever 

ijv.tr — ijv.tr'- In such case, since both sequences satisfy the bounds of Proposition 16.31 we conclude 

||P* - p:,\\^, \D* -DU< Cr cN . (6.34) 

Now, given <r, a' and n there exists N > n and a, a' s.t. <xj = crj for all i < n, similarly for the primes, and 
ijv,* = ijv,5-'- Indeed, just take N = n+\m(n, o J )—m(n, o~)\ and tr- = 0, a[ — 1 for i > n in case m(n, a')—m(n, a) > 
and vice versa in the opposite case. Hence, we conclude that (|6.34p holds for all a, a' and TV and the claim follows. □ 

Since the claims of the Propositions 16.21 and 15731 are uniform in the initial time to £ [to,2to] as well we have 
a- independent limits P* to = P t *,D* to = D\ Q . Let us also denote the to-dependence of the times (|6.33[) explicitly 
as i n , CT ,to- Let U c [tq, 2tq] be the set of to s.t. there exist a, a 1 and n, n' such that i n , CT ,fo = tn',o-',T - By the same 
reasoning as in the proof above, one deduces that on the set U, the limits P t * , D^ are constant. 
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Lemma 6.2. //log 2/ log f is irrational, then U is dense in [tq,2tq\. 

Proof. Consider \ogt n ,a,t — —2 log A + log to + 2nlog7? + 2m log 2. If log 2/ log t is irrational, then the map x — > 
x + log 2/ log^ mod 1 has dense orbits and the claim follows easily from this. □ 

Using a density of orbits argument as in the previous Lemma it is easy to see that there exists T < oo such 
that the set of times £n,er,t with to G U is dense in [T, oo). Along any sequence of times from this dense set, the 
limit (|2.29p exists and is independent of the chosen sequence. By the strong continuity of the Z t (which will be 
easily derived in Lemma 110.11) . the function (J2.29I) is continuous in t and hence the limit is independent of ^,to- 
This completes the proof of Theorems [T] and [2] Theorem [3] follows analogously by considering sequences of times in 

(i2~m 

7 Flow of T 

As announced in Section lB.1.31 we prove the Induction step for T. The following convention applies to Sections [5] 
and HO We always assume that the induction hypotheses Propositions 16.11 and 16.21 are satisfied on scale n (we then 
prove them on scale n + 1), and we do not repeat this assumption in the statements of all lemmata and propositions. 
Throughout this Section, as well as in Sections [8] and [9l we use the conventions on £, eo, e^o that were explained at 
the beginning of Section 16.11 

7.1 Powers of T n 

We show how to bound powers of T n and we abbreviate T — T n whenever no confusion is possible. This bound will 
be an important ingredient of both induction steps for T n and G c n A . 

Lemma 7.1. T — T n satisfies 

' dxe 10 ^^\\(S e T m )(x)\\ w <C, l<m<£ 2 (7.1) 

supe* > ' nlx[ \\(StT*)(x)\\v<C (7.2) 

X 

20-7Q I I 

supeW 1 * 1 \\T m (x)y <Cm- d / 2 , 20 < V™ < I (7.3) 



Proof. Let us first suppose max(20, -^) < y/m < I. By analyticity 



e^ lxl T m (x)= fdpf(p + i^^e x ) m e ipx (7.4) 

./ \/m 



where e x = x/\x\. To evaluate this integral let Xn(p) be the indicator function for |Rep| < p„. On its support (and 
with |Imp| < 70 we have f (p) m = e mf ^R(p) + ((1 - R(p))f(p)) m and using (|6T5| 

d PX n(p)\\f(p + i^e x ry < Ce^o [dpe~-i D " m P 2 + (±b n ) m p d n < C(m" d / 2 + «5 m ) (7.5) 

with 8 < 1 (since p n increases sublinearly with n ). 

To perform the integral with 1 — % n , we use first (I6.13|) and m > 20, and then m > ^, to obtain 

dp(l-x„(p))||T" Tl (p + i^e :c )||^ < b"Vol(T„) < C m r* mn £ dn < b^ 2 (7.6) 
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where Vol(T„) = (2ir) d £ dn is the volume of the d-dimensional torus T n = £ n T d . Combining these bounds we 
we get a bound 0(m~ d / 2 ) for (17. 4[) which implies (|7.3[) for max(20, i|p) < y/m. Taking m — £ 2 and using that 
(S e T rn )(x) = £ d T m (£x), we get (TT2l) . 

To get (|7.1[) . we first note that, from (|7.3p we get immediately 



i22a| x | 



dxeV^ w \\T m (x)y <C 



(7.7) 



!()(/ 



After a change of variables x = £y and using i^=j£ > IO70 (j7.1[) follows for ^/m > max(20, ^). For 20 < y^ < 
we get ||T m (x)||^e^ 7o|a:| < C from eq. (|JTT4l . This settles ([731) and, upon scaling, also flTT]). □ 

For later purposes we register the following consequence of (|7.5j) and ([7. 6ft : 

IK^Tf )(z) - f„ +1 (x)|^ < cr cn e - 20 ^l-l (7.8) 

where (p := p + 20i7oea;) 

f n+1 (x) ^ e- 2 °-">^ [ dpR n (p/£)e e2 ^^e^ x Xn(p/£)). (7.9) 



7.2 Contribution to T n+ \ from G C A 

We recall the expression for T n+ \ in terms of quantities at scale n; 



T, 



n+l 



M 



£ T T >St 

Ae<s(i T ,) 



® G>U ® r„(r) 

Ae.4 T^Supp^t 



(7.10) 



and below we define £"„ as the second term on the RHS. 

We now derive the necessary bounds on E. This bound is the only place in the induction step T — > T" in which 
we need information on the cumulants G C A . 



Lemma 7.2. Let 

Then 

and hence in particular ||£Vi|| 7o < Ce n 



E n :— T n+ i — Se[T n ]. 
svLpe 10 -<^\\E n (x)y<Ce n 



(7.11) 
(7.12) 



Proof. We (again) abbreviate T = T n , E — E n and G A = G c n+1 A , We take £ > 10 and apply Lemma [5.31 to eq. 



%> 



® G A <g> 

AeA r^Supp^t 



T(r) 



< n ii^n 7 „nii TiJ| 



70 



(7.13) 



1070 AeA 



where the product Y[j runs over all discrete intervals J in the sets J T '\Supp^l. That is; decompose the set / T '\Supp„4 
into a union of sets J of consecutive numbers ("discrete intervals") such that no two of them are consecutive. By 
invoking Lemma l7.ll we bound 

||T |J| || 70 <C, since |J|<^ 2 (7.14) 

The number of stretches is at most 1 + |Supp.4| < 2|Supp.4| and hence we can bound 

\\E\\ Wj0 < E CISUPPAI II H G ^lli07o < C f^ 2 (7.15) 

Ae<B(I T ,) AeA £ " 

where the second inequality follows from Proposition 16.21 and the fact that YIaga l-^l ^ s a ^ ^ eas ^ 2- 
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To get a bound on sup,,, ||i?(a;)||^e 7 ' a: ', we must proceed differently; we split the contributions to E = Ej + Ejj 
in (|7.10p into those where there is at least one T, namely when Supp.4 ^ I r /, and those for which there are no T's, 
i.e., Supp.4 = I T i . For the second term, En, we use that 

l|£//||io 7n < £ C^^Y[\\G c A \\io, <C e2 l^/: (7.16) 

Ae < B(I T ,),S\ippA=I T , AcA 

where £ 2e originates as a crude bound for the number of partitions of £ 2 elements (bounding the sum over A G 
*B(/ T ')). Further, remark that En is translation invariant and use the a priori bound on the supremum norms given 
by 

sup\\Ei I (x)\\ w e )/ ^<£ dn [ dx\\E II (x)\\^ x ^£ dn \\E II \\^ (7.17) 

This bound originates from the fact that the 'unit cell volume' of the lattice X„ is £~ dn . Using now the bound on 
II-E77II1070 derived above, we get 

sup \E u (x' - x)\e 10 ^ x ' x '^ < £ dn (C£ 2 e n y 2 < Ce £- cn£2& < e„ (7.18) 

x,x' 

where the second and third inequality follows from e n — eo£~ an (by Proposition l6.2l) . 
We now deal with Ej. Here one uses the bound 

mpWE^-x)]]^^'-' 'I < J2 ]J\\G c A \\, ( S up\\T( x '- x )\\^ x - x '\\ (7.19) 

x AQ<S(I T ,),SnppAjtI^ AzA ^ x ' x ' ' 

< C*e„sup||T(a;)||^e 7o|:E| < C*e„ (7.20) 

X 

The second inequality is obtained by summing A as in (J7.15I) . To get the first inequality, first apply point 2) of 
Lemma 15.21 in all s-coordinates and then apply point 3) of the same Lemma in the x-coordinates. □ 



7.3 Analysis of T n+ \ 

We complete the induction step for T by establishing 

Lemma 7.3. Provosition \6.1\ holds on scale n 4- 1. 

Proof. Whenever confusion is excluded, we abbreviate T = T n , R — R n . From Lemma 17.21 we get 

f n+1 (p) = f n (p/£y 2 +E n (p), with sup \\E n (p)\\ 9 < Ce n . (7.21) 

Let first |Rep| > p n +i- To study the first term in (|7.21|) with Proposition 16 . 1 1 we need to consider two cases: 
1. |Rep| >£p n - By (jeTTS)) 

\\Tip/efU<l£ <~b„ +1 . (7.22) 

Now use Ce„ = Ce £~ &n < \b n+x by taking e < c(£) to get (f6~T3l) for n+ 1. 
2- pn+i < |Rep| < £ p n . By Proposition 16. II 2) we have in this region 

f(p/£f = R{p/£)e e2f "^) + (i _ R(p/l))f(p/if (7.23) 

To control the first term on the RHS, we use (ETLU]) . (|6~T3T> and (|6TT5|) to get 

\\R(p/£)e e2 f^W\y < e-^"P" +1 e c ^ 2 < cr *s(n+i) < h +1 _ (7.24) 
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The second term on the RHS of (|7.23|) is bounded by 

||(1 - R(p/t))f(p/tf\\ 9 < b£ < ib n+1 (7.25) 

Again, we may assume Ce n < -|b n +i and so (|6.13|) holds for n + 1. 

Let then |Rep| < p n +i- We apply spectral perturbation theory, see Lemma TA. II Referring to the notation in that 
lemma, we have A = f n+1 (p) and A = T n {p/£f , a = ^ 2f ^ p ' l \ P = Rn(p/f) and A x = B„(p). 
Let us first consider p in a wider strip |Rep| < pn+i and |Imp| < 270. In this region by (|6.15|) . 

|oo] - |e^»^)| > e-i^+i-^o > cTt^+i) (7.26) 

whereas from (I6.7P 

mo-a Po||^<bf = (|fr^f". (7.27) 

This is no bigger than |ao|/2 say and so (|A1[) . the first condition of the Lemma [A. 11 holds. Because of ||Ai|| < Ce n 
and (|6.6p . the condition (|A13|) holds if 

Ce„ < |oo| 

which is true for small eo ■ Lemma IA.1I then implies that the isolated eigenvalue persists and 

; /n+i(p) _ e ^ 2 /„(f) < Ceni ( 7 ,28) 

\\R n+x {p)-R n {p/l)\\^ < Ce n (\a \-\\Ao-a Po\\ w -Ce n )- 1 <Ce r n ' i £i^ n+ ^ (7.29) 

so the first claim in (|6.9|) follows. Furthermore, 

||(1 - R n+1 { P ))f n+1 { P )\\y < C\\R n+1 ( P ) - Rn(p/e)\y + Cb* + Ce n < b„+i. (7.30) 

It remains to iterate (|6.8[) . The analyticity of f n +i follows immediately. To derive the symmetry properties of /„+i, 
consider a lattice symmetry O. We start from the symmetry property (cfr. Section [57 



T n+1 = Ad(Vo 1 ))T„+iAd(T/ ) (7.31) 

By the result in Section 15.3.31 if follows then that 

f n+1 (0- l p) = I^ f n+1 {p)I p ,o (7-32) 

which implies that f n +i(p) = f n +i(0~ 1 p), and hence the Taylor expansion of f n +i(p) around p = contains no 
odd powers. 

Furthermore, by (17726)) and ([7728|) . we have f n +i(p) = ^ fn(jp/£) + j(j>) with 

sup |j(p)| < Ce„^"(" +1 ) < V^^ i5 " (7.33) 

|Imp|<27o 

The diffusion constant D n+1 is defined as the quadratic term in the Taylor expansion: D n+ i5ij = ^-V»Vj/ n +i(p = 
0). By Cauchy, 

|^„+i-A 1 |<2 7 o" 2 sup |j(p)| (7.34) 

|p|<2 70 

and the second claim of (|6.D[) follows by the bound (|7.33[) . Furthermore, for |Imp| < 70 we have 

\f n+1 (p) + D n+lP 2 \ < \t 2 f n (p/£) + D nP 2 \ + \j(p) - fJ^PiPjCViVjMO)] (7.35) 

i,j 

< (£-l£-f » + 67 -3 sup |j(g)|)|p| 3 (7.36) 

|Imp|<27o 
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which proves (|6.8j) at n + 1 because of (|7.33j) . 

Eq. (I6.11[) follows since T n necessarily preserves the trace of density matrices and Trp = J2 s es p(p = 0> s )- 
The bound (I6.12p iterates due to (|7.29|) and will be established for n = in Section fl 1.31 



We complete the proof of items 2) and 3) of Proposition 16.11 by noting that the constants C in the bounds (|6.6l 
16.121) may be chosen uniform in n by invoking (|6.9[) for all n' < n. To get item 1) of Proposition 16. II . note similarly 
that the proofs of 2) and 3) give for all p with |Imp| < 70 

\\f n+1 (p)-f n ( P )\\^<cr cn 

with c > 0. Since this holds with n replaced by n' for all n! < n, T ra+1 (j>) can be bounded by an n-independent 
constant. The bound on T„ +1 (0,p) propagates trivially due to T n+1 (0,p) = T n (0,p/£). 

It then remains to prove item 4) . For this we use the L°°-bound on E(x) (Lemma 17.21) . and the result (|7.8[) . 
Together they show 

e 10 ^ x \\\T n+1 {x) -f n+1 (x)y < cr cn . 

Consider now f n+1 (x) defined in (|73|) . From (J6T51) and ([7T55|l we get \Re£ 2 f n (p/l) + D n p 2 \ < C on the support of 
Xn- Hence 

e^ x \\\f n+1 {x)\y<CD- d l 2 sup \\Rn(p)h- 

|Rep|<p rl ,|Imp|<7o 

These estimates show 4) holds uniformly in n. 

This finishes the proof of the induction step n — >• n + 1 for Proposition 16.11 on the condition of the bounds of 
Proposition 16.21 at level n. □ 



8 Flow of G C A : Linear part 



.4* 

We will now analyze the recursion relation for the correlation functions G A given in p. 341) . In the present section 
and in Section [SJ we drop the subscript n on operators, writing e.g. G C A = G c n A ,T — T n and we denote operators 

on scale n + 1 by a prime, e.g. G A = G° +1 A 

We split the recursion relation into a linear and a nonlinear part (although these adjectives do not fit literally). 
Consider the terms in the sum in (|3.34p with A a collection that consists of just one set, namely A — {A}. We use 
the notation A — > A' to indicate that A contributes to G A , in this sense. In other words, 

A -> A' & AeI A ' and Vr' e A' : (A n I T >) ^ 0- (8.1) 

The linear RG flow is the contribution of such A, with the additional restriction that \A'\ = \A\: 



A:A-*A',\A'\=\A\ 



G C A ® <g) T{t) 
reI A ,\A 



J.2) 



The aim of this section is to state good bounds on \\G A , H n || lo-yo; such that the induction hypothesis Proposition ^. 21 
can be carried from scale to scale if G A is replaced by G A lin . We call the remaining terms in the sum (J3.34I) nonlinear 
and we treat them in Section |H1 By inspecting Proposition 16.21 it is easy to understand the qualitative difference 
between the linear and nonlinear contributions. Recall that in Proposition 16.21 we claim a bound proportional to 
e« . The nonlinear contributions to G A , carry at least \A'\ + 1 powers of e n and this gives us some leeway. However, 
the linear contribution has just \A'\ factors of e n and this forces us to do a careful analysis to be able to get a bound 

\A'\ 

with e„ +1 . There are two mechanisms at work that help us here: 

• Rescaling of time: Since the microscopic times described by the macrotimes r on scale n + 1 are £ 2 larger 
than those on scale n, the blow-up factor dist(A) Q for A — > A is roughly speaking £ 2a (\ A I -1 ) times larger 
than dist(A')". 
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• Ward Identities: The Ward identities from unitarity and (in the case 0\ = 02) reversibility allow to improve 
our estimate on each term in the sum on the RHS of (|8.2[) . The gain from each of the Ward identities is 
roughly £ _1 . (this will be explained later) 

The main force working against us is 

• Entropy factors: For each A', there are £ 2 \ A ' sets A such that A — > A' , \A\ ~ \A'\. However, due to the 
summability over A with min A fixed in Proposition 16. 2[ only one factor I 2 (corresponding to min A) , shows 
up in our estimates. 

The gain from the Ward Identities is explained in Section l8~T1 Then, in Sections 18.21 and l 8.3[ we prove Proposi- 
tions 18.11 18. 2\ 18.31 which we present now. 

Proposition 18.11 provides a bound on G A , lin which is valid regardless of whether p^ is an equilibrium state or 
not. For this reason we call this bound 'non-equilibrium'. In what follows, we denote by \x\ the largest integer 
smaller than x. 

Proposition 8.1 (Non-equilibrium linear RG). 

]T ^t{A l r\\G%^A\ 10l0 <C\ k J l - ia r^ k ^{\og c (.) 1 k>2 (8.3) 

A' CN:|yl'|=fc,min A' = 1 

J2 dist^riiG^aiio^^-^iog^) (8.4) 

A'CK:\A'\=2,minA' = l 

Note that the bound in the case \A'\ = 2 is only useful when a > 1/2, since for a < 1/2 the power of t on the 
RHS becomes nonnegative. The next two propositions are only meaningful in the equilibrium case 0\ = 0% = 0. 
Proposition 18.21 deals with correlation functions G c n A with A 3 0. 

Proposition 8.2 (Initial time linear RG). Assume 0\ — 02 — 0, then 

]T dist(A') tt ||G^ ain || 1070 <C fc 6^„r 2 «LVJ(iog c £), k >2 (8.5) 

A'CN :|A'| = fe,minA'=0 

£ dist(Ari|G^ lin || 107o < Cele hn r^ 2a \\og c l) (8.6) 

A'Cl ! i Q :\A'\=2,minA'=0 

The improved bound in Proposition 18.21 (compared to Proposition 18- 1[) is due to the fact that the macroscopic 
time does not get rescaled, and hence there is essentially no entropy factor from sets A containing 0. 

Finally, we state a bound that applies to bulk sets A, but that is better than Proposition 18. II since it allows to 
treat the \A'\ — 2 correlation function for a < 1/2. 

Proposition 8.3 (Equilibrium linear RG). Assume 0\ = 02 — and that < a < 1/2. Then 

J2 dist(Ar\\G c A , x J\ 10l0 <Ce 2 n r 2a (log c e) (8.7) 

A'CN:|A'|=2,minA' = l 



The improved bound in Proposition 18.31 is due to the fact that we can use two Ward identities instead of one. 

If one neglects the contribution of the nonlinear part in the recursion relation (|3.34[) . then the above bounds 
establish the induction hypothesis Proposition 16.21 Indeed, for 1/2 < a < 1 (the case 0\ ^ 02), Proposition 18.11 
implies the bound (|6.17j) of Proposition l6~2l provided that 

l-2a < -2a (8.8) 

(l-4a)-2a[(fc-3)/2j < -ka, k>2 (8.9) 
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and this is indeed satisfied for our choice a = a/2 — 1/4. For 1/4 < a < 1/2 (the case fix = fa), we need (combining 
Propositions 18.11 and I8.3|) 

-2a < -2a (8.10) 

(l-4a)-2aL(fc-3)/2j < -ka, k>2 (8.11) 



to satisfy the bound (|6.17|) and (by Proposition I8.2[) 

-2a < -2a -ai (8.12) 

-2a[(k-l)/2\ < -ka-ai, k>2 (8.13) 

to satisfy the bound (J6.18I) . One inspects that a = a/4 — 1/8 and a\ = 1/4 do the job. 

8.1 Contraction from the Ward Identities 

Consider the connected correlation function G A £ Si 'a (we have dropped the subscript n again). As explained in 
Section SHI it satisfies the Ward identity 

dx' T (P®l...®l)G c A {x' A ,XA) =0, r = maxA (8.14) 

where P is any projector in & of the form P = \fj)(ls \ with X) s es M s ) = 1> as introduced in Section |4~T1 We will 
use this identity to derive a bound on the operator 

(T m ® 1 ® . . . ® 1) G^ (8.15) 

where T m acts on the leg of the tensor product indexed by r — max A. The relevance of (18.151) is that it comes up 
after writing out the RHS of (|8.2[) explicitly, i.e. performing the contraction Ta> ■ However, for the purpose of the 
present section, one can just accept (|8.15p as the basic object of study. The reason why one can get a bound on 
(|8.15[) is that, formally speaking, T m has a right eigenvector with eigenvalue 1 given by ls = l[ s es ] (in particular, 
it is independent of a;) and it contracts vectors orthogonal to l,s to size ~ to -1 / 2 . The condition (|8.14[) ensures that 
G c n A as a function of x' T is orthogonal to ls and hence (18.15[) should decay ~ mT x l 2 . This intuition is captured by 

Lemma 8.1. Let K £ Si be a kernel satisfying the 'Ward Identity' 

fdx' PK(x',x) = 0, (8.16) 

Then, there is an operator V m ,i & Si satisfying 

\\S e V mno/2 \\ W7o < C1 ° g( i- +m) , for 1 < m < t\ (8.17) 



and such that 

|T m in<|T/ mi7 ||Xe^'-*l|, (8.18) 

where we recall that \L\ is the operator whose kernel is the absolute value of the kernel of L, and the inequality in 
the last formula is between kernels. 

The proof of this lemma is postponed to Section 18.1.11 

Now it is straightforward to bound (|8. 15|) (after rescaling), indeed, let 7 = IO70, then 

\\Si(T m (g>l(g>...(g>l)G c A \\ W70 = \\\(T m (g>l(g>...(g> 

< (\\V m , lo/2 <E>l<Z>... 

^ C ^ C t +m) \\G< 



1)<-M 1070/1? 




(8.19) 


® 1) (G c A e^ lx ^ 


-*-*'"* Al )\\\l0 10 /i 


(8.20) 


|| iOTO i TO 




(8.21) 
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The first inequality uses (|8.18[) with the coordinates ZMmaxAt z A\m&xA kept fixed. The second inequality uses (|8.17|) 
and the definition of the norm || • || 7 . 

Let us now try an analogous trick based on the Ward identity (|4.31[) which we rewrite here with the help of the 
projectors P? := \^){l So \ and P ref = |// of )(l 5o |; 

dx a G\ lT} {x a ', x' b ; x a ,x b )(l ® P?) = f dx a L T {x a \ 4; x a ,x b )(l ® P rcf ) (8.22) 

It will be used to get smallness for the right action by T m i.e. for G?j -, (1 ® T m ) by a similar reasoning to the 
one given above for the Ward identity from unitarity. The operator T formally has a unique left eigenvector with 
eigenvalue 1, given by fj,r„ as defined in 16. Ill and T m contracts vectors orthogonal to /j,t„ to size ~ to -1 / 2 . Hence, 
if (|8.22|) would hold with the RHS replaced by and P$ replaced by R(0), then wc could repeat the reasoning of 
the unitarity Ward identity. However, the projection P 13 — \fi^){ls \ is ^/eo-close to P(0) = |Mt„)(1s d I' see &3M in 
the induction hypothesis. Actually ^t„ converges to /i" asn-) oo, but this is not even needed here. The presence 
of the RHS in (|8.22|) introduces an additional error term which is small because L T consists of correlation functions 
that include r = and therefore contract faster. 

Lemma 8.2. Let K,K £ ffl be kernels satisfying the 'Ward Identity' 

J dx K(x', x)P 13 = J dx K(x', x)P rci (8.23) 

Then, there are operators W mn ,W mn £ S% satisfying 

\\StW m , yo/2 \\ 1(ho < Cm-^ 2 logm + C\\R(0)-P^\y, (8.24) 

II^W m , 7o/2 ||ioy < C (8.25) 

for 1 < m < I 2 , and 

\KT m \ < \Ke^ x - x '\ \W mrl + \K<f\ x - x '\ \W mrl (8.26) 

We will use the Lemma to derive a bound on the correlation function (T m + <Ei 1)0^ , (1 (g) T m -). Its kernel is 
bounded 

| (T m + ® 1) G\ lT} (l(g>T m -)\< V m+tl \e^<- x ^G c {ltT} (1 ® T m -)\ (8.27) 

< V m+n (\e^\<- x -\ + \<- x ^G c {XtT} \W m _ n + \e^<-^L T \W m _ n ) (8.28) 

where we used first (|8.18|) and then (|8.26|) . Collecting the bounds on V, W, W, we get 

\\St(T m +®i)G c {1<T} (i®T m -)\\ Wj0 
< c^|||e^K--iq 1)T} (i0T™-)|| 1O7o/ , 

< C l -^± (wGl^W^^i 1 -^^ + ||P(0) P% s ) + ||L||«J) (8.29) 

The operator norm ||£ T ||7 can be bounded by norms of correlation functions: 

\\L T \\-j < l|G{ 01/r} || 7 + ||G{ 0;T }|| 7 + ||G{ 0t _ 1} || 7 (8.30) 

by (|4.3ip and Lemma \5. 31 (1). The upshot of the calculation in (|8.29[) is that all terms between the large brackets 
on the last line are smaller than \\G% -,\\ wy (which would be the resulting bound if we had used only the unitarity 

Ward identity). As anticipated before Lemma l8~2l the smallness comes either from the factor og _ff^ , the difference 



||i?(0) — P \\<g an d the correlation functions contributing to L T that are small because they involve the initial time 
and therefore contract faster, see Proposition [821 
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8.1.1 Proof of Lemma 18.11 

We start from (|8.16[) . i.e. J dx PK(x,xo) — for any xq. This implies trivially that 

T m K(x\x Q ) = fdxe-^ x ' xal (T m (x' ,x) - T m (x' \x )P) K(x,x Q )e^ x ~ x ° l (8.31) 

We now choose the projector P to equal R(0) (the spectral projector of T(p = 0)) and we define V mn € 3fc 

V m ,-y(x', x) := sup e - T|x - Xo1 {T m {x - x') - T m {x - x')R(0)) (8.32) 

x 

where the absolute value | • | is applied to a kernel in s, s', that is, for an operator D in Sf we set \D\ to be the 
operator with kernel |D|(s', s) := \D(s' , s)\. Note that V mn is translation invariant by the translation invariance of 
T m . From ([831]) and (|Q2]l . we get the inequality (|8J8]) of Lemma O 

It remains to establish the bound on V mn , i.e. (|8.1T|) . It suffices to consider y/m large enough only, we take 
yfm > 80. Since the operator V ma is translation invariant, we can use Lemma 15.41 and hence it suffices to show 
that 



Indeed, we then get 



supewW||y miTD/2 (0 )a; )|k < Cm-VW+Vlogm. (8.33) 

10-fO 20-10 - 



\V m ^o/2\\w~, /i < s m I dxe y 



1)1x1 < Cm- 1 / 2 lo 



gm 



(8.34) 



which implies (|8.17[) by scaling. 

To prove (|8.33[) we split V mn (x',x) = 1^ 7 (i',i) + V r ' J ' l , y (x',x) where the terms are defined by restricting the 
supremum in (I8.32[) to ^-\x — x \ < logm and ^|x — x \ > logm respectively. We bound 



^"-(i 



^ W IC 7o/2 (0,z)l < e^ 1 * 1 sup e-^\'-">\\T m (x)-T m (x )R(0)\ 



x -to\x-xo\>iogm 



1 2°ai x | 

m 



-"-(i 



T m {x)\ +sup-^e'^ lx ' x ol e v^ |xo1 |T m (x o )i?(0)| 



xq| and i/m > 80 in the last step. The desired bound (|8.33[) 



where we used the triangle inequality \x\ < |xo| + 
(without the log) now follows from f|T.3[) . 

Next, we consider V' m , 2 . The argument for showing (I8.33|) for V' m , 2 is analogous to the proof of Lemma ITT 
and we will use some notation from that proof. We start from the momentum space representation 



T m (x',x)-T m (x',x o )i?(0) 



dpe 



ip(x— X ) 



f m ( P ) - f m (p)R(0)e- ip( - x - Xo) 



.35) 



and proceed as in (I7.5[) and (I7.6I) to conclude that up to terms exponentially small in m (and n) it suffices to bound 

f dpe mf " (p) e ipx (R(p) - R{0)e- ip( - x - Xo) ) X n(p) (8.36) 

JT„ eg 

where p = p + —pS-e x . This is bounded by 

C f dpe-^ mp2 (\\R(p) - R(0)y + |1 - e-W x — >|)*n(p). (8.37) 

using the bound (|6.15l) . Since 

|1 - e-^-^l < C{\p\+"f /V^~) logm 

the second term on the RHS of (J8.37I) contributes 0(m~( d+1 ^ 2 logm). For the first term use analyticity of R(p) in 
a ball of radius 70 at the origin to get 

\\R(p) - fl(0)||* < C(\ P \ + 7 o/V^+ l|Rc P |>c 70 ) 

and hence the first term on the RHS of (|8.37[) contributes (D(m~( d+1 ^ 2 ), as well. 
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8.1.2 Proof of Lemma f872l 

We start from the Ward Identity 

/ dx'K(x", x')pP = / dx'K(x", x')P lcl (8.38) 
Then 

dx'K{x",x')T m {x',x) (8.39) 

dx'K(x", x')e^ x "- x '\ (T m (x', x) - P^T m (x", x)) e^ 1 *""*' 1 (8.40) 

+ / dx'K(x",x')e~< lx "- x ' l P Icl T m {x",x)e-^ x "- x ' 1 (8.41) 



We now define the operators W m ,-y, W m , 7 by specifying their reduced kernels 

W m ^{x',x) := sup|T m (x',x)-P' 3 T m (a;",a;)|e- 7|2: "- a; ' 1 (8.42) 

x" 

W mn {x',x) := snp\P Icf T m {x'\x)\ e -'' lx "- x ' 1 (8.43) 

x" 

(the absolute values on the RHS are again meant as absolute values of kernels on A„ x A„ , see the remark below 
(|8.32p ) Note that both operators are translation invariant. By completely analogous reasoning as the one leading 
to (prfe]) . we get the inequality (f8T26|) . 

Now to the bounds on W mn , W mn - If P^ is replaced by R(0) in the definition of W mtl (x' , x), we get the first 
term of the bound by an analogous proof as that of Lemma 187X1 The error term due to R(0) — P^ is estimated by 
C\\R(0) — P^\\y where the constant C originates from bounding SiT m with the help of Lemma l7.1l Similarly, the 
bound on W m . 7o /2 is immediate from Lemma |7. II 

8.2 Preliminaries for the induction step G C A — > G C A , 

In this section, we gather some tools for the proofs of Propositions 18.11 18.21 18.31 We state the main bound that 
we will use to control the sum over correlation functions at the lower scale. We abbreviate k := \A\ = \A'\ and 
m + :— max Ia> — max A; 



l|G'A',lml!l07 < 2 , 

A->A',\A\=k 



SiTa' 



Ga ® T(t) 

TEI A ,\A 



1070 

k 

< £ \\(T m +®l®...®l)G c A \\ Wlo/i \{\\T^-^- l \\ Wlo/i , 

A-+A',\A\=k j=l 



, fc log(l + m+) 

A^A',\A\=k 



£ cfc! f=ll^l^ ( 8 - 44 ) 



On the second line, we set the dummy tq = or, if t\ = 0, then tq = —1. The bound is obtained by proceeding as 
in the proof of Lemma [7721 bounding powers T m by Lemma l7~T| except the power T m + , which is bounded by using 
Lemma 18. 1[ as explained in the beginning of Section 18.11 

Next, we introduce some useful notation. We let h(A) stand for a function of the ordered times Ti, . . . , Tk of A 
with the property that, for min^4 > 0, 

h(A) = h(A + t) (time-translation invariance) (8.45) 
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and satisfying the normalization conditions, uniformly in Tfc,Ti, respectively, 

J2 H A ) < 1. and Yl h ( A ) - l ( 8 ' 46 ) 

0<Ti<T 2 <...<T k :T k fixed Ti<T 2 <...<T fc :Ti fixed 

Note that the first condition follows from the second by time-translation invariance. In particular, we will use that 
the function of k — 1 variables obtained by performing the sum ^ r . r <T h(A) or J2 T -t >t M^) satisfies the 
same conditions, and hence we will also call it h(-). This is how h(-) will enter: we write 

l|Gl lin ||io 7o < £ ^l log (j_t^) dl st(A)-"fe(^)- (8-47) 

by using (|8.44p . Now our task is to estimate the RHS, which is done by elementary analysis in the next sections 
(note indeed that all operators have disappeared from the RHS). Below, we consistently use the notation ti,T2, . . . 
for the ordered elements of A and t{ , 7% , . . . for those of A' . The variables t\ , Ti , . . . range over the sets I T > ,I T > , . . . 

8.2.1 Change of variables 

We need to sum (|8.47j) . Hence we will have to convert the factor dist(A') Q into dist(A) Q . In doing this, we will gain 
small factors of £~ 2a . Indeed, we find the inequality 

dist(Ti^) a <r 2a dist(A) a , A^{Ti^},T^-r[>l (8.48) 

That is, if the macroscopic times are not neighbors, then we gain a small factor. However, even if the macroscopic 
times are neighbors, then we still get a small factor from every second difference of macroscopic times, i.e. 

dist(A') Q < Cl A 'l (r 2a ) [1 ^ r ~ Li dist(A) Q , for any A -> A' with \A'\ = \A\ > 3 (8.49) 

Indeed, for all 1 < j < \A'\, either \tj — Tj +i \ > ||rj — rj +1 | or \tj — Tj_i| > ||rj — rj_ x |. We will refer to these 
bounds as 'change of variables'. 

8.3 Proof of Proposition 18.11 

To prove Proposition ^. 11 we need to control the sum 

]T dist(AT £ 1 di S t(A)-*h(A) (8.50) 

A'llA'l^fe^axA'fixed A-^A' V T fc 

which emergetQ from the LHS of (|8.3|8.4|) by (|8.47|) . We will do this by using the change of variables and the 
integrability of h(-), both introduced above. Note that we dropped the factor log (1 + m + ) since it can always be 
estimated by C" log I which suffices for our purposes. We also did not write e n since this factor just goes through 
all estimates. 

8.3.1 The case \A'\ = 2 

Let us assume first that t[ < t' 2 — 1, then (J8.50I) reduces to 

6 In fact, Proposition 18.11 demands that we keep rain A 1 fixed. However, by time-translation invariance, the claim 
with maxA' fixed is equivalent 
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E ( 1 + T 2-r[) a J2 MTi,r 2 )(l + r 2 -r 1 )-«(max/ r , -T2 + 1)- 1 / 2 (8.51) 

< C£' 2a E E Mn,r 2 )(max/ T , -rs + l)- 1 / 2 (8.52) 

< cr 2a J2 {™&xir>-T2 + i)- 1/2 = cr 2a Y,x~ 1/2 <ct 1 ~ 2a (8-53) 

T 2 ei, x=\ 

T 2 

To get the first inequality, we used the change of variables formula (|8.49[) , the second inequality follows from the 
properties of h(-), i.e., from ([8.46J1 . 

Next, let us treat the case r^ = t[ 4- 1, then (|8.50j) gives (we can bound dist(A') Q by a constant here) 

C E /i(n,T 2 )(l + r 2 -r 1 )^(max/ r ,-r 2 + l)- 1 / 2 (8.54) 

rie/ T j,r 2 e/ T , 

< CJ2{1 + r 2 - min/^)- Q (max/^ - T 2 y x l 2 = ^ x~ a (l + ^ 2 - 20~ 1/2 < Cf- 2a (8.55) 

In the first inequality, we used (1 + r 2 — ri)~ a < (1 + t 2 — min/ r <) _Q , and we bounded the sum over t\ by (|8.46|1 . 
The last inequality follows easily after splitting the sum in 1 < x < £ 2 /2 and £ 2 /2 < x < I 2 . Putting the two 
contributions together, we get C£ x ~ 2a as a bound for (|8.50[) in the case \A\ = \A'\ = 2. 



8.3.2 The case \A'\ = 3 

We start with the subcase t{ = t' 2 — 1 = t 3 — 2 which turns out to be the most tricky one. First, we dominate 

E dist(T 1 ,T 2 ,T 3 )- tt (l + max/ T , -ra)- 1 / 2 ^^!,^,^) (8.56) 

< E dist(r 2 ,T 3 )- Q (l + T 2 -min/ r ,)-«(l + max/ T ., -r 3 )- 1 / 2 Mr 2 ,T3) (8.57) 

T2.3EA-' 

T 2,3 

by (1 + t 2 — ri) _Q < (1 + t 2 — min/ r /)~" and the properties of h, i.e. ()8.46[) . Then, we change variables 

x — r 2 — min/ T j, y = max/ r ' — t 3 , < x,y < £ (8.58) 

such that we have to bound 

e 2 -i 

J2 F(x, y), F(x, y) = h(x, 2£ 2 - y - l)(x + iy a {y + iy l ' 2 (2£ 2 -x-y- \y a (8.59) 

x^y— 

where we used the invariance of h(-, •) under joint translations of its arguments. 
To perform these sums, we define the sets 

C a n ■- {(x,y)\(l~6 nfi )2 n <x<2 n +\ (1 - 5 n , )2 n < y < 2/U 2 } (8.60) 

C b n ■= {(x,y)\(l~6 nfi )2 n <y<2 n +\ (1 - S nfi )2 n < x < 2/U 2 } (8.61) 

C c := {{x iy )\2£ 2 -x-y<2/3£ 2 } (8.62) 
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And we will split 



53 n*, y) < E E F & y) + E E F (*> y) + E F (*> ») ( 8 - 63 ) 

x,y=0 n=0 C a n=0 C b C c 



sup i 



for n = 0, . . . n* with n* the smallest natural number such that 2™ > £ 2 . To perform the sum in C%, we bound 

((2£ 2 -x-y- l)- a (x + l)- a {y + 1)~ 1/2 ) = cr 2a (2 n y 1/2 ~ a , (8.64) 

we use properties (|8.46|) to perform the y-sum and we bound the s-sum by 2™ , the 'width' of its domain. This yields 

J] F(x, y) < Cr 2a {2 n ) 1 - 1 ' 2 - a (8.65) 

The sum in C h n is done analogously, except that now the i-sum is controlled by (j8.46[) and the y-sum by its width 
2". On C c , we can bound (x + l)~ a (y + l) -1 / 2 < £~ 1_2q , and then the sum is done straightforwardly as 

J2 F ( x >y) < ^ 2a E Hx,2l 2 - y)(2£ 2 - x - y - iy a (8.66) 

C c %,y=0 

e 2 -i 
< C x - 2a J2 (t 2 - x)- a < l 1 -^ (8.67) 

x=0 

The sum over n — 1, . . . , n* yields 

e 2 -i 

E F fav) < i 1 ~ ia + r 2a J2(' 2n ) 1 ~ 1/2 ~ a ( 8 - 68 ) 

x,y=0 n>l 

< Cr 2a {2 1/2 - a ) n " < Cf~ 4a (8.69) 

Let us now look at the case where one pair of t[,t 2 ,t^, is consecutive. If 1 + t[ < r^ = 73 — 1, then we get 
a factor t~ 2a from (1 + t 2 — ti) q by change of variables, we integrate t\ by using the normalization of h(-) and 
we can treat the remaining times r^Tg as outlined in the case \A'\ = 2, gaining an extra factor £}~ 2a . The other 
possible case with one consecutive pair, i.e. 1 + t[ = t' 2 < 73 — 1, can be related to the previous case by symmetry 
considerations, using the translation invariance of h(-). Finally, the case where no pair is consecutive is of course 
analogous to the corresponding case with \A'\ = 2: One gets £~ 4a from change of variables, and the remaining sum 
over 71,2.3 yields I by using the normalization of h(-) and the factor (1 + max/ r < — t^)^ 1 ' 2 . 

Hence in all subcases with \A'\ — 3, we obtain a factor i 1 ^ 4 ". 

8.3.3 The case \A'\ = k>3 

We perform the sum over 71,72, ••■ ,Tk-3- As argued in (|8.49[) , we get at least one small factor £~ 2a from every 

second sum. This yields the bound £~ 2q L 2 J . The last three sums are then bounded by C£ 1 ~ 4a by repeating the 
analysis of the case \A'\ = 3. 

Combining all cases, we get the claim of Proposition 18. II 

8.4 Initial time linear RG: Proof of Proposition 18.21 

We start as in Section 18.31 from the bound (|8.50p . Now we perform the sum over times n, . . . ,Tk starting from Tk, 
keeping in mind that t± = 0. 
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In the case k = \A'\ > 2, we get a factor £~ 2q Lt-J by the change of variables, and the sum is performed by 
using the summability of the function h(ri, . . . , Tfc) with t\ = fixed. In contrast to the proofs above, this suffices 
since t\ is indeed fixed here: t\ — 0. For the case k = \A'\ = 2, we perform the sum over t 2 by splitting it in the 
regions r 2 > £ 2 /2 and t 2 < £ 2 /2. For t 2 > £ 2 /2, we can use the change of variables, and hence we get 

(ESHD < c ^ r 2a h(o,r 2 ) + c J2 (^-^r 1/2 (i + T 2 )-^(o,T 2 ) 

t 2 >£ 2 /2 l<r 2 <e 2 /2 

< cr 2a + cr 1 (8.70) 

This finishes the proof. 

8.5 Equilibrium linear RG: Proof of Proposition 18.31 

We have to perform the sum 

£ dist(T{,T 2 n|Gf r ,^ lin || 107o (8.71) 

m+ = max/ T ' — T2, TJi- = T\ — min/ T ' (8.72) 



Let us abbreviate 



Then we obtain 

N/ oc' I, ^ V- ^ log(l + m+) ^ log(l + m_) 

ll G {^,^},lmlll07o < 2^ ° (1 + m + )V2 I (1 + m _)l/2 ^{n.^lU 

{ti,t 2 }-v{t{,t0 

" ) "ll ( ^{o,i,T 2 -ri+i}ll7o + l|G{ 0]T2 _ Tl+1} || 70 + ||G| 0/r2 _ Ti j|| 7D J (8.73) 

where we used the time-translation invariance property G c < T 1 = G c , x T _ T +1 i , the bounds f|8.29[) and (18.30[) . and 
the bound (1021) for \\P? - R{0)\\^. 

Next, we perform the sum over r{ of the RHS in (J8.73I) . This RHS is split in four terms. The r(-sum of the last 
three terms is bounded brutally by C£ e„ei )n by using the integrability of the function h(-) to perform the sum over 
r 2 — T\ and estimating the sum over t\ by £ 2 . By the smallness of ei jn , this bound is sufficient for Proposition 18.31 
Next, we focus on the ("r-J-sum of the) first term in (|8.73|) . It is of the form 

J2 dist(r{,r 2 ) Q J^ C(l + m + )- 1 / 2 (l + m_)- 1 / 2 dist(r 1 ,T 2 )- Q /i(T 1 ,r 2 ) (8.74) 

r l :T l< r 2 {t"1,T 2 }->{t^T^} 

To evaluate this sum, let us first consider the case t[ < t 2 — 1. We set 

z := t<z — n, x = max/ r ' — r 2 (8.75) 

and we estimate (|8.74|) restricted to t[ < t 2 — 1 by 

00 £ 2 -l 

Cr 2a ^/i(T 1 ,r 1 +z)^(l + a;)- 1/2 (l + (a: + z)mod£ 2 )- 1 / 2 < Cr 2a \og£ (8.76) 

z=e 2 x=o 

We used the change of variables to get the factor £r 2a and the Cauchy-Schwarz inequality to estimate the i-sum. 
Restricting (J8.74J) to t[ = t 2 — 1 and setting y := n — min/ T ' , we get 

e 2 -i 
C Y^ h(l,f-x- y)(l + x)- l ' 2 {\ + y)- 1/2 {£ 2 -x-y- \y 2a (8.77) 

a:,y— 

This sum is analogous to the one treated in Section [5 .3. 21 the only difference being that one exponent is 1/2 instead 
of a. The multiscale treatment can be copied without changes. The result is that the sum is bounded by Ci~ 2a 
and this yields Proposition 18.31 
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9 Flow of G C A : nonlinear part 



A 

In this section, no information on T n is needed, except for the bound on T™ from Lemma 17.11 Starting from the 
induction hypothesis on the cumulants at scale n, we deduce an estimate on the nonlinear part of the contribution 
to scale n + 1. We do not distinguish between bulk and boundary terms, and hence we have icNo throughout. 
We drop the scale subscript n and we mark operators on scale n + 1 by a prime, as explained at the beginning of 
Section |U 

As was explained in Section [S] as well, the nonlinear contribution to the RG flow is defined by excluding from the 
sum (|3.34p the contributions of A = {A} (a single set) with \A\ = \A'\. Hence we still need to study the remaining 
terms 



^A.nlin : ~ &£ 2-~t ^ A ' 



A e <3(Ia>), Qa' (-4) connected 
12a€a\ A \>\ A '\ 



<g> G C A ® <8> T(t) 

AeA r^Supp.4 



(9.1) 



that is, cither A has more than one element, or if it consists of one element A, then \A\ > \A'\, and this is combined 
in the condition J^AeA 1^1 > l-^'l- ^ ur resu lt is 

Proposition 9.1. Fix an exponent a with d < a and a constant c v > 0. If c v is chosen small enough, then 

sup J2 (e')" |A,| dist(A') Q ||GL llin || 107o <l, with e' = c~ l r & e, (9.2) 

T ' A'Cfin-A'Br' 

We remind the conventions introduced at the beginning of Section |6~T1 £ _1 and eo should be chosen sufficiently 
small compared to constants like c v . Proposition 19.11 establishes (the nonlinear parts of) Proposition 16.21 Indeed, 
for a > 1/2, it suffices to fix a < a < a whereas for the case a < 1/2, we observe that we can choose a such that 
additionally ka + a.\ < k"a holds for k > 2. 

To derive Proposition 19. 11 one can ignore the Ward Identities completely since they can at best give a factor £ 2 
in each term of the sum, and such factors are irrelevant in view of the fact that we have extra factors of e due to 
the condition J^AeA l-^-l ^ l-^-'l ( cn ~- ^ ne discussion at the beginning of Section[8]). 

Consequently, we will use the following basic bound for the cumulants. In fact, a special case of this bound 
appeared already in the proof of Lemma 17.21 

Lemma 9.1 (A priori recursion relation). 

l|Gl||io 7o < E eC ' SupMI II H^IU (9-3) 

A&B{I a i),9 a i (A) connected AeA 

Proof. We take £ > 10 and apply Lemma [5731 to cq. (|3.34p : 

% 



® G C A ® T(t) 

AeA -r^Supp.A 



< n ii G Aii7 nii TiJ| ii7o (9.4) 

1070 AeA J 



where the product Y[j runs over all discrete intervals J in the sets I T > \ Supp^4, cfr. the proof of Lemma 17.21 By 
invoking Lemma r7.ll we bound 

||T |J| || 70 <C, since|J|<^ 2 (9.5) 

The number of discrete intervals J is at most 1 + |Supp.4| < 2|Supp„4| and this yields the claim. □ 

9.1 Sum over connected coverings 

Our strategy for bounding the sum over polymers A in eq. (|9.1|) consists of the following splitting. For each A G A, 
we let S(A) C A' be the macroscopic domain of A, that is 

S(A) = {T'e A',AnI T >^$} (9.6) 
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Note that any collection A of sets A induces a collection S = S(A) with elements S(A). We call a collection of 
sets connected whenever it can not be split into two collections whose members are mutually disjoint. For any 
A € Q5(/a')> the connectedness of the graph Ga'(A) implies that S = S(A) is connected and that Supp<S = A' . We 
call the set of connected collections £. 

With this terminology, (|9.3[) can be written as 



\G 



A'IU07o 



* E E II ll<^M° |A| 

5G£,Supp5=A' yle»(/ A / ):S(7l)=5 Ae.4 



(9.7) 



< 



e n* 

5££,Supp5=A'SeS 



E H^II 70 e C|A| 



(9.8) 



where the function F\{x) := X)^°=i xP = x /(l — *) appears because there can be more than one set A E A such that 
A ->S. 

Let us now determine how this bound can be modified if we restrict the sum in (|9.3[) to those entering in the 
nonlinear RG, i.e. in (|9.1|) . If |«S| > 2, then the condition J2agA l-^-l > l-^'l * s a l wavs verified. If |<S| = 1, i.e. there 
is one S = A', then there are either at least two A £ A such that A — > A', or there is only one but it satisfies 
\A\ > \A'\. Hence we get 



Gc 



A',nlinl|107o 



< 



E II Fi [ E \\ g a\\ ,.' 

S6£,Supp5=A'|5|>lSe5 XA^-S 



.C\A\ 



+ 2^1 H^AII '■> 

A^A> ,\A\>\A>\ 



c .ll e C \ A \ 



F 2 ( E II? 

\A->A' 



All7o e 



(9.9) 



(9.10) 



where F 2 (x) =J27= 



p=2 ' 



The following lemma shows how to control sums over A — > S that will appear in the evaluation of the above 
expression 



Lemma 9.2. 1) 



2) Abbreviate 



Then 



E e^ldist^niG^IU < (Ce)^ 

A^S 

v n (S) := (r & ep s \dist(Sr E H G All7oe C|A| . 

A-> 

SUp E MS) < CI 2 



A^S 



f)2+2a 



SBt' 



3) 



sup E (^ & e)~l s ldist(S) tt E ll G All7o eC|j4 ' < <^ 2+3a e. 

T ' 59r' A-vS,|A|>|,S| 



Proof. To get (J9.11I) . we recall the function h(-) introduced in (|8.47[) . and we bound the LHS by 

E( Ge )"'E E h(A)<f(CefK ri=mmS 

k>\S\ nG7 T / A:\A\=k,minA=T 1 

using the summability of h, i.e. (|8.46j) and bounding the sum over t\ = minj4 by £ 2 . 



(9.11) 

(9.12) 
(9.13) 

(9.14) 
(9.15) 
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Now to (|9.13j) : we rewrite 

Y, U "(S)< E (^ & ty lSiA)l £ lAl dist(S(A)) a {dist(A))- a h(A) (9.16) 

S3t' A:A(M t ,^$ 

Assume there is a r' G S(A),t' ^ min S(A), max S(A), such that the set A n / T < contains only one element. Then 
we get a factor £~ 2a from change of variables, and if there are p such r', then we get at least [(p + 1)/2J of these 
factors. From elementary considerations, 



2([( P +l)/2\)>p>\S(A)\-(\A\-\S(A)\)-2 (9.17) 



Hence, we bound 



S3t' A:AnI T ,^ 

< £ 2a E (e£ a ) kA - ks £ (& - a)ks E h ( A ) 

k A >k s >2 A:AnI T ,^H,\A\ = k A 

< £ 2a+2 V fc j4 (e^) fcA -' cs f( <i - a )' cs < C£ 2+2a (9.18) 

fcA>fc s >2 

The second inequality follows by setting k$ = \S(A)\, Ua = \A\. To get the third inequality, we bound the sum over 
h(A) with A 3 t by \A\ times the sum over h(A) with min A fixed, and then we proceed as in the proof of, f|9. llj) . 
This yields (|9.13j) . To obtain (|9.14|1 . one repeats the calculation with the only difference that the sum over kA,ks 
in (|9.18j) is restricted to fc^ — 1 > k$ > 2. 

D 

9.2 Proof of Proposition 19.11 

We perform the sum 

E (e'r m dist(Ar l|Gl,, nlin ||io 7o < I + II + HI (9.19) 

A'3t' 

where the three terms on the RHS refer to the three terms in (|9.10[) . The second term can be bounded immediately 
with the help of (|9T4|) (with S = A'), yielding 

II < C£ 2+3a e (9.20) 

We estimate the term III, with x := J2a-^A> WGaWio^^^ ^ (Ce)^' 1 by (J9.11I) in Lemma [9~2l In particular, we 
have x < 1 — c, and hence F 2 (x) < e 2 Cx with F 2 as in (|9.10p . since \A'\ > 2. Therefore, 

III < Ce 2 J2 (£T lA ' l dist(A') a J2 l|G^|! 70 e c|A| < e 2 £ 2+2a (9.21) 

A'3t' A—>-A' 

where we used (|9.13|) in Lemma I9T21 (recognizing the definition of v n (S) with S = A'). Now we turn to 

I< E (n- lA ' l dist(AT £ II E \\G A h e ClAl (9-22) 

a'bt' see,SuppS=A' ,\s\>i ses A-ts 

where we dropped the function i*i(-) at the cost of increasing the constant C, using, as in the treatment of term III 
above, that its argument is smaller than 1 — c. First, for any S € €. with SuppS = A' , Lemma fB.ll implies 

dist(yl') Q < [7 dist(5) Q (9.23) 

ses 
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such that ()9.22j) is dominated by 

P C\A\ 



sec,\s\>i,S\ippS3T> ses \a^s 



(9.24) 



By substituting (J97L2"T) for each S e S in (|9~24| . and setting N(S) = J2ses \ S \> 

(15^21 < J2 (g/)Af(5)-|SuppS| ffl (2+2a)|5| ~Q c |g| r (2+2a)^ (g) (g25) 

SeC,|S|>l,SuppS3T' "^ 'ses 

By the connectedness of S, we have N(S) — |Supp<S| > |<S| — 1, and combined with |<S| > 2, this yields r(S) < 1. 
Hence we are left with the task of estimating (the constraint |<S| > 1 can now be dropped) 

E II ^k-^^v n {S). (9.26) 

see,SuppS3-r' ses 

The tool to perform this sum is the bound ^2g 3T i £~& +2a 'v n {S) < C from Lemma W7Z\ Indeed, sums over collections 
of connected sets, as in (J9.26I) . can be handled conveniently with cluster expansions, with the bound from Lemma 
19.21 playing the role of the Kotecky-Preiss criterion. We state the relevant cluster expansion result in Proposition 
IB. II in Appendix [El A basic corollary, eq. (|B5j) . implies that (J9.26J) is bounded by a constant that can be made 
smaller than 1 by choosing c v small enough. 

Collecting the bounds on the three terms I, II, III, Proposition 19. II follows 

10 Estimates on the first scale: Excitations 

In this section, we prove bounds on the correlation functions G A , i.e. the claims of Induction hypothesis 16.21 for 
n = 0. We treat the bulk correlation functions (0 3 A) in Section 110.21 and the boundary correlation functions 
(0 G A) in Section 110.31 We derive an explicit representation for the operators G C A that is based on the Dyson 
expansion. This representation will also be useful to analyze T for n — in Section QTJ Since we start so to say from 
scratch, the first part of our discussion is in finite volume A and we perform the thermodynamic limit in Section 
110.2.11 

10.1 Dyson Expansion 

We recall the reduced dynamics of the system at microscopic time t 

Z tPs = Tr E [e- HH (ps ® f% f y tH ] (10.1) 

and the dynamics on the n = scale, T„ = o (introduced in Section [3]), it is related to Z t through To = Z\-2 io . In 
the next sections, we will develop an expansion for Z t and we will relate this expansion to the correlation functions 
G A . 

10.1.1 Derivation of the expansion 

We introduce the time-evolved interaction Hamiltonian 

fl>(«) : = e isffE E (<f>(q){W ® e iqX ® a* tq ) + h.c.) e~ isHE (10.2) 

q£A*,i=l,2 

and define the Liouvillians Lj(s) = ad(-H/(s)), is = ad(iJs) and Lb = ad(i/g). Let us also use the shorthand 
U s = e~ lsLs . Then the Duhamel formula 

e itL Ee -itL = Ut+ f dsU t ^s(-i\L I (s))c isLE e- 
Jo 
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yields upon iteration the Dyson series for the reduced dynamics Z t : 

Z t ps = X)(-A a ) m / dt 1 ...dt 2m Tr E [u t - t2m L I (t 2m )...U t2 - tl L I (t 1 )U tl (p s ®p% i )\ (10.3) 

m ^° 0<ti<-..<t 3m <t 

where the invariance TiE(e~ ltLE A) = Tte A was used. For m = 0, the RHS is understood as Utps- 
Next, we write the integrand using the formalism developed in Section [3. II 



Tr E 



U t - t2m Li(t 2m ) . . . U ta -tMti)U tl (pB <8> Pe ) (10-4) 



= T E [U t - t2m ®s Li(t 2m ) ®s • • • <8>s U t2 - tl ®s i/(*i) <8>s CtJ Ps- (10.5) 

We recall that the expectation E acts on an element of J?® <E)S?e and the contraction T : M® — > S?. As in Section 
13. 1[ it will be convenient to label the spaces Si in S?® by the times that occur in the corresponding operators. 
Given a 2m-tuple of times {t%, . . . ,t2m} = < C [0,£] let io = 0,t2m+i = t and denote the family of intervals 
Ji = [ti,ti + i],i — 0,2m by J{t). We will then index the space S?, where Lj(ti) lies by Mt t and the space M where 
U ti+1 - U lies by St Jt . Thus 

2m 

E[Ut-t 2m <S>sLi(t 2m )®s---®sUt 2 -t 1 <E)sL I (t 1 )® s Ut 1 }e ® S? u <E> S?j (10.6) 

and as before the operator T contracts the operators in the obvious chronological order, i.e. such that those in S?[o itl ] 
are on the right, then those in S?t 1 , then those in &H lt t 2 ], etc. 

Let {u,v} be a pair of (distinct) times with the convention that u < v. Then we define 

K u . v = -A 2 E [£,(«)) ® s Lj(ti))] , K u . v eS?<S)Sg (10.7) 

and we view this operator as an element in S? v <g> S? u . We also abbreviate U s i_ s by Uj with J = [s, s'] and we view 
Uj as an element in S?j. Since Lj is linear in the creation and annihilation operators, the Wick theorem yields 

^[Ut-t 2m ®sLi(t 2m )(g> s ...<S) S Ut 2 -t 1 <Z>sL I {t 1 )(g>sUt 1 }= V <8> K u>v ® C7> (10.8) 

where the sum on the RHS runs over pairings it, i.e. partitions of the times t\, . . . , t-zm m m pairs (u, v) with the 
convention that u < v. 

By plugging (|10.8[) into (|10.3[) . we obtain 

Z t = J2 / dti...dt 2 m J2 T'[^{u,v}^Ku,v^j&j{t)Uj\ (10.9) 

m ^°0<ii<...<i 2m <t 7rePairings(t) 

As before, we equip products of S%, e.g. as in (|10.6[) . with the norm || ■ || . 

10.1.2 A formalism for the combinatorics 

The integral over ordered t, together with the sum over pairings, tt, on the set of times, is represented as a combined 
integral and sum over ordered pairs (itj, Vi) with m, Vi £ M. + and i = 1, . . . , m, such that 

m < Vi, m < . . . < u m (10.10) 

This is done as follows. For any pair (r, s) £ tt, we let m = t r ,Vi = t s where the index i = 1, . . . , m is chosen such 
that the u* are ordered ui < it2 ■ • ■ < u m . We represent one pair (it*, vi) by the symbol Wj and the m-tuple of them 
by w. We call J7j the set of w such that Ui,v% £ J (for arbitrary m), and we use the shorthand 

dw:=^2 / d - / d ^ 1 [ti ( <«i]l[u 1 <...<ii m ] (10.11) 

Qj m ^° J- J'" 
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where the RHS is set to 1 for m = 0, corresponding to w = in the LHS. In what follows, we will often consider 
the ordered times t,u,v to be implicitly defined by w. For example, we will write J(w) instead of i7(|). The Dyson 
expansion in terms of the sets of pairs w reads 



Z t = 



dwT 



lw&wK w <g> U.j 
JeJ(w) 



(10.12) 



«[0,t] 



10.1.3 Connected correlations and the Dyson series 

To relate the previous sections to the setup in Section^ we need to discretize time and express the operators Ga, G a 
in terms of the Dyson series. 

Recall that N is the set of macroscopic times. To each macroscopic time, we now associate a domain of 
microscopic times, 

Dom(r) = [A- 2 t (r - 1), A^t] (10.13) 

To a set AcNof macroscopic times, we then associate the domain 



Dom(A) = (J Dom(r) 



(10.14) 



t£A 



We take t = X 2 toN. Then, a set of pairs w G £l[o.t] determines a graph Q(w) on {1, . . . , N} by the following 
prescription: the vertices t < t 1 are connected by an edge if and only if there is a pair w = (u, v) in w such that 



u G Dom(r) and v G Dom(r ) 



(10.15) 



(Note that there may be several such pairs). We write S\xpp(Q(w)) for the set of non-isolated vertices of G(w), i.e. 
the vertices that have at least one connection to another vertex. If w G ^DomA, than Supp(tJ(w)) is obviously a 
subset of A. In that case we write Qa{w) for the induced graph with vertex set A. The graphs G(w) in the Dyson 
expansion with support A give rise to the correlation function Ga of Section[3l and connected graphs Qa{w) give rise 
to the connected correlation functions. This goes as follows. Recall the collection J(w) of intervals determined by 
the times in w. Given a macroscopic time r G N we define J T {w) as the family of intervals {Dom(r) n J \ J G J7(w)} 
and set 

J A (w) := (J J T {w) (10.16) 

t£A 

Using the group property of U s and the definition of contraction we have 

U.j = i(U. h ®U. h ) 



for J — Ji U J2 with Ji, J2 consecutive intervals. Applying this to the intervals J in (|10.12j) intersecting more than 
one Dom(r) we get 



T 



® K w ® Uj 

wEw JeJ(w) 



T 



® K w (g) Uj 

wew J£J{i,n](w) 



--■ TV[o,n](w). 



(10.17) 



where we abbreviated [1, N] — {1, 2, . . . , N}. The tensor product defining V[i^](w) factors across the macroscopic 
times, i.e: 



V[i,N](w)= ® V Dom(r )(u;)(g)Vsupp(a(H))(^) 

T^Supp(g(w)) 



which can also be written as 



T^sup P (e(uj)) * 



(10.18) 



(10.19) 



where Qi are the connected components of Q. 
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As in Section (|3.5j) . we may perform the time ordered contraction T in two steps, first within the time intervals 
Dom(r) and then contracting the rest: 

(10.20) 



N 

T = T <8> T T 

T = l 



where 



i:tiSDom(r) JeJ T (w) 



The beautiful formula (|10.20[) is not a typo but a consequence of the fact that we defined T both as contracting 
microscopic times and intervals, and macroscopic intervals. Writing as before Ta = ®t£A%-i e q- (jlO-ISp leads then 
to the expansion (|3.6p with 



G A = 



d-W l[Supp(G(w))=A]T~A 



® K w <g> Uj 

wew Jej A (w) 



'Dom(A) 



and 



dw T 



® K w ® Uj 

wGw JGj T (w) 



'Dom(-r) 



(note that r on the RHS is arbitrary) whereas (|10.19j) gives 



a 



d«Z \q a (to) connected] IA 



® K w ® Uj 

w£w j£j A (w) 



'Dom(A) 



(10.21) 



(10.22) 



(10.23) 



10.2 Bounds on bulk correlation functions 



In this Section, we state and prove Lemma 110.11 which is actually the induction hypothesis 16.21 for bulk sets A (not 
containing 0) and on scale n = 0. Boundary correlation functions (A 3 0) will be treated in Section [10.31 



10.2.1 Bounds on the Dyson expansion 

As a first step, let us derive a term by term bound on the Dyson expansion and perform the thermodynamic limit. 
Let us write the operator K UjV in (|10.7[) explicitly. Using the field operators ( I2.15[) we can rewrite (|10.2[) as 



Hi(s) = ^3 w ® 1* <8> ®(x, s) 



(10.24) 



Let us use the notation ad(A) = A — A 1 where ^4° is left multiplication and A 1 is right multiplication (by A). 
Then, (fTOJ)) becomes 



k u , v = -x 2 J2 Y, (-i) a+b ( ab (x-y,v-u)(w®t x ) a ®(w®t v ) b 

x,V a,be{0,l} 



(10.25) 

(10.26) 

(10.27) 

(we used translation invariance in time and 0(d) invariance in x). 

Eq. (|2.17p (thermodynamic limit for Q implies the kernel of K UjV has a pointwise limit as A /*■ Z . For the 
remainder of the present section, we use the notation K u ^ v and h(s) (introduced below) both for A finite and 
A = Z rf , indicating differences whenever necessary. Since the kernel of 1" is diagonal in the coordinates x, v (recall 



where 



By (|2~T6| 



C b (x -y,v-u) = Tr E ($(x, v) a *(y, ufp^) 



C( x ,t) = c w (x,t) = c 11 ^,-*) = r(x,-t) = c(x,t) 
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that z = (x, v, 7], ei,, en) £ A$ and note that this v has nothing to do with the time-coordinates u, v used below), we 
get 

\\K U ^\\ 7 < CX 2 (10.28) 

for all 7. In fact, using the time decay of C m assumption fA"! and denoting 



we get h(s) < C for A finite and 



X 2 h(v-u) :=\\K U , V \\ 



ds(l + \s\) a h(s) <C, 



2070 



for A 



(10.29) 



(10.30) 



Lemma 10.1. The sums and integrals on the RHS of (|10.12[) , p0.21[) p0.22[) and (|10.23[) converge absolutely. For 
example, the series defining Z t is bounded by 



dw 



^[0,t] 



T 



® K w ® Uj 

w£w JeJ(w) 



< 



20 70 



3 CX ' A finite 
s cx2t A = 7L d 



(10.31) 



In particular ; the limits of (|10.12[) , (I10.21[) (j!0.22[) and p0.23[) as A /* Z d exist. Moreover Z t is strongly continuous 
in t. 



Proof. The following reasoning applies for finite and infinite A alike. Using Lemma 4.3. we get 



T 



<S> K w ® Uj 

wEw JeJ(w) 



^ n ii^ii-r x n \\ u j 



7 w£w Jej(w) 

and using standard propagation bounds for the lattice Laplacian A, we have 

||t/,|| 207 o<^ 2 ^l 

(recall that we treat 70 as a constant). Hence the LHS of (|10.31[) can be bounded by 



f dm. II Ce x ' c \ J \x J] X 2 h{v~u) 

q J£J[o,t](nL) ™ 6 ™ 

e x2ct J2 f du(f[C J dvMvi-Ui)) 

ox't sr ( x2ct ) 



E 



' 0<Ml<...<M m <t 



II) 



dsh(s) < 







cx ' A finite 



^cx-'t 



A = Z d 



(10.32) 

(10.33) 

(10.34) 
(10.35) 
(10.36) 



To get the second last inequality, we first performed the t^-integrals, and then we estimated the w^-integrals by the 
volume of an m-dimensional simplex. In the last one we used eq. (|10.30l) and h(s) > C for finite A. For any finite 
t, the term-by-term convergence of the series follows by the Vitali convergence theorem, since each term converges 
pointwise and the series is summable uniformly in A. For the series in (|10. 211 [TO. 221110. 2"3)l . similar reasoning applies. 
Strong continuity of t i— > Z t follows from strong continuity of s 1— > U g . □ 



The following estimate is an immediate consequence of (|10.23|) and bounds as in the proof of Lemma 110.11 

Lemma 10.2. 

l|G^||io 7o , <e c|j41 I dw\\X 2 Ch{v~u) (10.37) 



Q 



Dom(A) 



w£'W_ 



w spans A minimally 

The statement 'w spans A minimally" means that Qa{w) * s connected and that no pair can be dropped from w 
without losing this property. In particular, this implies that Ga{w) is a spanning tree on A. 
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Proof. By (|10.23|) and Lemma TlO-H we have 

||Gl||io 7o < J dw H Ce A2c|J| x H X 2 h(v - u) (10.38) 

^Dom(A) 

Sa (™) connected 



v, J£Ja(w) w ^m. 



and, since |J7a(w)| < |ty|+|A|, we may dominate the integrand by e '" 4 ' E[«;ew A 2 C/i(?;— w), since |Dom(A)| = A~ 2 |^4|. 
Next, we state an appealing estimate was the main motivation for encoding the pairings -k in the pair-sets w. For 
any (integrable) function / on ftr>om(A)j we have 

I dw\f(w)\ < f dw' f dv/'\f(w!Uv/')\ (10.39) 

S Dom(A) ^Dom(A) £do„i(A) 

Qa (w) connected w' spans A minimally 

To realize why this holds true, choose a spanning tree £? for the connected graph Qa{w) and then pick a minimal 
subset w' of the pairs in w such that Ga(w!) — S? ■ Since, in general, this can be done in a nonunique way, the 
integrals on the RHS contain the same w more than once, and the inequality is strict unless / is concentrated on 
minimally spanning w. 

We apply this inequality with f(w) := YiwGw ^ 2 Ch(v — u) to (|10.38[) (with the integrand dominated as indicated 
there). The resulting integral over dw/' can then be performed similarly to (|10.36[) . yielding e c ' >. This proves the 
claim. D 

We now derive the main result of the present section 
Proposition 10.1. For A small enough and with €q — C|A| a , 

J2 eo lAl dist(Ar\\G c A \\ 107o <l (10.40) 

A C N : min A = 1 

Proof. Each w that spans A minimally determines a spanning tree on A. Hence we can reorganize the bound (|10.3T[) 
by first integrating all w that determine the same tree. This amounts to integrate, for each edge of the tree, all 
pairs (u, v) that determine this edge. Furthermore, we use that 

dist(A) Q < Yl ( x + l T ' - r l) Q ; for an y spanning tree ,9 on A (10.41) 

{r,r'}e£(.T) 

where £ (,5?) is the set of edges of the spanning tree £? (see Appendix [Bl for a simple proof). We arrive at the bound 

dist(A) a ||G^|| 107b < e c ^ J2 II ^( r ' T ') ( 10 - 42 ) 

span, trees ,5* on A (t.t')gS (S) 

where, for r < r' ', 

e Q (r,r') :=A 2 (l + |r-r'|) a / du f dvh{v-u), (10.43) 

J Doni(r) J Dom(r') 



and e a (V,r) :— e a (r, r'). Next, we establish 



E Ur,r')<C\X\ 2a , (10.44) 

t'SN\{t} 
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by using the bound (|10.30j) on h from Lemma Tl 0.1) 

e Q (r,r + l) < X 2 C f ° ds sh(s) < \ 2 C(2\- 2 to) l - a < C\\\ 2a 
Jo 

A - rto Co 

J2 ea(r,r') < A 2 C(t |Ar 2 r" f du f dv(v - u) a h(v - u) < C\X\ 2a 

T ' >T+1 A- 2 (r-l)io u+X-^lo 

Starting from the inequality (|10.42[) . we bound the sum in (I10.40|) as 

J2 e^ Al dist(Ar\\G%\\ 1070 < J2 (C\\\ 2a )-\ £ ^ J] e a (r,r') (10.45) 

A CN trees ST on N: (r,r')e£(^) 

minA = l Supp(^) 3 1, \&\ > 2 

where we denote by Supp(^) the vertices that have at least one edge and £{3T) is the set of edges. This sum can 
be controlled relying on the bound (|10.44[) . This is a special case of Lemma \B. II in Appendix [Bl where (110.44[) plays 
the role of the Kotecky-Preiss criterion (|B3I) . k = 1, and the edges {r, r'} of the tree 2T play the role of the sets S. 

n 
10.3 Bounds on boundary correlation functions 

In this section, we work in the equilibrium case /?i = 02 — /3, since only in that case the boundary correlation 
functions are relevant. We recall from Section 14.2.21 

Z tPs = Tr E [e- itH D ld (ps <E> p^)D; d e itH ] (10.46) 

which differs from the reduced evolution Zt by the fact that we included the Radon- Nikodym derivative D Td . 
For a dense set of p G S%\ (Jf) , we can write a formal Duhamel expansion 

CO p 

D ldP D:. d = e AF ^Y, I dp 1 ...dp m \ m L I (/3 m )...L I (p 1 )p (10.47) 

m=0 0</3i<...</9 m </3 

where the term with to = is defined to be 1 and 

L I {p i )p = - l -[e^+^) Hie -^+^) ip]+ (10.48) 

with [B, A\ + = BA + AB. We will now combine this expansion with the Dyson expansion for the unitary evolution 
to obtain an expansion for Z t . For this it is convenient to set in eq. (|10.47[) fa =U — with ti G [— fa 0] and 

Li{t) :=Li{t + p) for te [-0,0]. 

We also generalize 

Uj :=C/./n[o,t] for Jc [-/Mi- 
ami 

K u>v := (il[„> ] + l[ u <o])(il[«>o] + l[»<o])E[i/(v) ®s L/(u)] (10.49) 

This bizarre looking formula simply takes care of the fact that the Dyson expansion for the Radon Nikodym derivative 
does not have factors of i. In particular, the above definitions of Uj,K u ^ v reduces to the ones given previously in 
Section HO. 1.11 in the case J G M + ,0 < u < v. For a set of pairs w G Q[-0 t i\, we let now J{w) be the induced 
collection of intervals partitioning the interval [— (3, t] instead of the interval [0, i\. 
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It is now straightforward to check the formal expansion 

(10.50) 



e AF ^ I dwT 



<8) K w <E) Uj 



The factor e AF ^- ) can be determined from the relation 

1 = TiiZ pf] = Tr[Z (l [;c=0] z/ cf )] (10.51) 

where p l g f and v rci were defined in Section [4.2.11 (note that l[ x =o\v is a function on Ao) and the last equality 
exploits the fact that p™ ! and v rof do not depend on x e Xo and also the operator Zq is translation-invariant. The 
advantage of the rightmost term in (|10.51[) is that Zq acts on an operator in £%i(J?s) that is strictly localized on 
the lattice and in particular it has a limit as A /* oo. Let us first define 

h(s) = sup 1 1 K u<v \\ 2 o 10 (10.52) 

— @<u<v:v — u— s 

Then, Assumption |A1 still implieaj that J R ds(l + s a )h(s) < C. To see this, let us inspect how the operator K u/U 
written out explicitly in (|10.25|) . gets modified for negative u,v. 

1) Depending on a, & € {0, 1}, there are minus signs that do not affect our bounds. 

2) The summand in (|2.16p gets multiplied with e ±rulj ( q \ < r < f3. Since the function ujj is smooth and bounded, 
this does not affect the decay properties of K U;V for large v — u. 

3) The operators 1% in (|10.25[) should be replaced by e rHs t^.e~ rHs , < r < (3/2. By a simple propagation estimate 
in imaginary time, cfr. (|10.33|) . the || ■ ||207 -norm of such expressions is uniformly bounded by a constant. 

By the above remarks, it now follows that Lemma 110.11 still holds if we replace h —¥ h. In particular, we can 
bound the expansion in (|10.50[) as in Lemma 110. 1[ we control the expansion of the rightmost term in (|10.51[) and 
hence we obtain the thermodynamic limit of the number e ^' and the equilibrium state r/^, and we prove that the 
operator Ti(0) defined in Section [4.2.2l is bounded and has a thermodynamic limit, too. Finally, from the expansion 
we get as well that 

sup fe 20 ^l t 'l|/(s) - // of (s)|) < CX 2 (10.53) 

ses ^ ' 

with (J,p(s), // ef (s) defined in Section 14.2.1 1 and we recall that s — (v, r], eL,en). 

Let us next generalize the expression for the correlation functions G C A . For the macroscopic time r = 0, we 
define Dom(r) := [— f3, 0]. Then the equalities (I10.21[) and (|10.23l) remain true without any changes and Lemma 
110.21 remains true upon replacing h —¥ h. 

Finally, we prove the bound (|6.18[) (Induction hypothesis I6.2J) at scale n = 0. 

Proposition 10.2. Recall that eo = C|A|" and ei : o = C|A| 2 ~ 2 ". For A small enough; 

Y, e^ |A| dist(^) a ||G^|| 207o <e I ,o (10.54) 

AcN 
min A = 

Proof. The proof mimics that of Proposition 110.11 We first derive 

dist(A) Q ||G^|| <e c l A l J2 II S «( T ' T ') ( 10 - 55 ) 

span, trees 3? on A (t,t')(z£(S') 



7 Of course, this is true provided that we treat j3 as a fixed constant and we do not attempt to get bounds that 
are uniform as f3 — > oo. 
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where e Q (r, r') was defined in (|10.43[) and this definition carries over to the case where one of r, r' equals 0, provided 
that we replace again h by h. However, the bound (|10.44|) can now be improved as 

__ f C|A| 2q r^O 

E e a (ry)<\ (10.56) 

the improvement for r = is due to the fact that the length of Dom(r) is f3 < C for r = instead of toA~ 2 for 
r^O. In the sum over trees (|10.45[) . there is now always one edge of the form {0,t} for which we have the bound 
CX 2 = CegCijO instead of Ce 2 ,, and this accounts trivially for the improved bound on the RHS of (|10.54[) . □ 

11 Estimates on the first scale: reduced evolution T 

The analysis of the operator T = T n -o proceeds in three steps. Note that T depends on the coupling constant A both 
through the strength of the interaction and through the definition of the time scale A _2 to, since T = Zx-^u- In a 
first step, accomplished in Section ril.il we prove that, as A — > 0, the operator Z\-2 t can be replaced by the Markov 
semigroup e*(~ lA L s+ M ) with M a dissipative operator. This is the 'Markov approximation' that we referred to 
already in Section 12.4.31 This Markov approximation is well-known in the literature since the work of Davies [7] 
and we describe it mainly for reasons of completeness. 

In the second step, in Section HT.21 we argue that the Markov semigroup has good properties, corresponding more 
or less to the requirements that Prop [B~3l imposes on T. The reasoning in that section is standard in the analysis 
of (classical) linear Boltzmann equations, involving standard tools as Weyl's theorem and the Perron-Frobcnius 
theorem. Finally, in Section 111.31 we argue that T inherits these properties from the Markov semigroup. This 
follows immediately by perturbation theory 

11.1 The weak-coupling limit 
11.1.1 The generator M 

We first define the generator M and we investigate its properties. Let U® stand for the propagator XJ t with A = 0, 
that is; 

jjo = e -itL. Pto) Lsp . n = ad (# spin ) (n.i) 

and we will again write Uj to denote Ug,_ s for an interval J = [s, s']. We define 

\ 2 G° V _ U := T [K a , v ® U^ v] ] (11.2) 

A 2 G„_ U := T[K UtV ®U [uM ] (11.3) 

Since U® is a finite dimensional unitary tensored with the identity, acting only on the spin degrees of freedom, 
ll^lli < C uniformly in £,7, and (|10.30[) implies 

/•OO 

/ d*(i + Mr ||G°|| 2070 < a (11.4) 

Jo 

Next, let P e be the projector to the the eigenspace corresponding to the eigenvalue e E £ = er(L sp i n ), cfr. eq. (|2.21[) . 
Then the connection between the operators G® and the generator M is via 

M = V/ dsc- [sE P e G° s P e (11.5) 

see Jo 

Thus we also have 

||M||20 70 < C. (11.6) 
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11.1.2 Emergence of M from the Dyson expansion 

We recall the expansion (I10.12|) for the reduced evolution Z t : 



Z, 



dwT 



[0,t] 



® K w ® Uj 

wEw J£j(w) 



We will show that the leading contribution to Z t in the above integral comes from of pairs 

fi[0,t] := i^ G °[0,t] : V i < U i+1 * = 1; • • • ) N - !}• 

Set ilh t , = £l[o.t] \ ^rn t] anc ^ fix 7 = 2O70 in Lemmata II 1.11 and II 1.21 First we have 
Lemma 11.1. 



(11.7) 



(11. 



dw 



[0,t] 



T 



i wew K w <E) Uj 

J£j(w) 



CVtl \\2ct 



<Ce UAt \X 



Proof. We proceed as in the proof of Lemma 110.21 and wc bound the LHS of (jll.9[) by 

e x2ct f i dw [] X 2 h(v - u) 



(11.9) 



(11.10) 



[o,t] wGw 



Every w £ flh t , has to contain at least two pairs (u, v), (u',v') such that u < u' but vl < v. Choose the first two 
such pairs and integrate over the coordinates of all other pairs, using that \\h\\i — J dsh(s) < 00 and proceeding as 
in the proof of Lemma 110.11 This yields the bound 



fTrToD < c A 2 *(c+Nli) gm 



(11.11) 



with 



q(t)= / h(v - u)h{v' - v!) 

and q(t) can be easily estimated by first performing the integral over v', which gives a factor \\h\\i, and then the 
one over u' , such that one gets (recall that a < 1) 

(11.12) 



(11.13) 



q(t) < \\h\\i\ 4 I dudvh(v - u)\v - u\ 

0<u<v<t 

< ||/i||iA 4 / du( max Is] 1 '"] [ ds\s\ a h(s) 

Jo \0<s<t-u J J 

< mJ fd S \ s rh( S ))(\hr\\\ 2a 



which yields the bound (|11.9p upon invoking Assumption |A"1 



(11.14) 



D 



Next, we organize the contributions from leading sets of pairs. We call Z® the integral over them, i.e. in eq. 
i)l 1 .T|) we replace ^[o,t] by f2? t , . Note that these can be written explicitly as 



00 „ 

Zj^ / A I Cu/ClU . . . U U3 ~ V2 Lj V2 — U2 U U2 — Vl \-t Vi —u\ ^ u\ 



(11.15) 



where G s and U s were defined above. We prove 
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Lemma 11.2. 

\\Z° - C <-- iLs+x2M ^\\ 7 < Ce A2ct |A| 2Q (11.16) 

Proof. Using integrability of G t and the bound (|10.33|) together with the product structure in eq. I|11.15|) , we get 
that the Laplace transform 

/>00 

Z° := / &te- tz Z G t (11.17) 

Jo 

is analytic in the half plane Re z > CX 2 and given explicitly by 

4° = (z + ilspin + iA 2 m- x ad(A) - A 2 G Z ) ' (11.18) 

where G z = L dte~ tz Gt and we recalled that L$ = L sp [ n + X 2 m~ 1 ad(A) . Using (|10.30[) we get G z is analytic in 
Rez > C|A| 2 too and bounded there by 

||G,|| 7 <C. (11.19) 

The same bound holds for ad(A). Recall that the spectrum of L sp i n is real and given by the set £ (|2.2ip . Hence 

||Z?|L < - , —r (11.20) 

11 zllf ~ dist(z,i£) v ; 

if Rez > C|A| 2 . By the resolvent identity 

Z° -(z + ii^)- 1 = A 2 4°(im- 1 ad(A) - G z )(z + LL^)- 1 (11.21) 

which implies, using (|11.20|l and Rez > C|A| 2 , 

\\Z z P s \^<C—±— (11.22) 

where we recall that P E is the spectral projector to eigenvalue e eigenspace of L sp i n - 
Inverse Laplace transform gives 

z o _ e (-ii s +A 2 M)t = f dze t^z° z -(z + iL s -X 2 M)- 1 ) (11.23) 

where C — {z : Rez = A|A| 2 } and A is taken large enough, A > G. By the resolvent identity we have 

Z° z -(z + iL s - X 2 M)- % = X 2 Z° Z (G Z - M)(z + iL s - A 2 Af)- J (11.24) 

and from (fPL~5j) 

M = J2^PeG° P 6 . 

The integrand in (J11.23J) is large when z is close to — ie, e € £ . Thus we localize the contour C — D s q£C £ such that 
on C £ , — ie is the nearest element of —i£. For z G C e we write 

G z -M = (G z - G z ) + (G° z - Gi) + (Gl - M) := D\ + Z>£ + D§. (11.25) 

and set 

X 2 Z°Df (z + iL s - A 2 Af y 1 := If(z) (11.26) 

Then, by ([TOT]) . 

Z° - e (-iL s +\ 2 M)t = <£ f dze tz (It(z) + II(z)+II(z)). (11.27) 
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|£>f|| 7 < / d^ft)|c- c 'l A l 2t -l| <sup ' n ■ ^ n ' ) / dth(t)(l + t) a <C\X\ 2a . (11.28) 



Recalling (J11.2JI and (|11.3|) and using \\U t - C/°|| 7 < e c|A|2 * - 1 together with Rez = A|A| 2 , A large enough, we 
bound 

'|e-Q|A| a «-lf 

where we used (|10.30p . 

To bound || -D§ ll-y we n °te that by (|11 .4[) G° is analytic in Rez > and by the same bound as in (|11.28[) . we get 
for z eC e 

/>oo 

||£>i|| 7 <C/ dth(t)\e~ zt -e [e \ <Cmin{l,|z + ie| 2Q } (11.29) 

Jo 

Since (|11.20|) also holds for the resolvent of — iLg + \ 2 M 

d*e'*(/f(,s) + /f(z))|| 7 <C|A| 2 / dze** -. (|A| 2a + minjl, |z + ie| 2 °})) < Ce CA2t |A| 2a . (11.30) 

c e Jc e l z + ie l 

Finally, to estimate If, from (jll.51) we have M = ^2 , P £ 'G® e ,P e ' and thus 

D% = Gl -M= Y. P ^ QPs. - E P c <%' P c 

Using (jll.22[) and the analogous estimate for P 6 (z + iLs — A 2 M) _1 we get 

||/f(z)|| 7 < CA 2 Yl (Iz + ieiWz + iEil)- 1 

and so 



|| / e tz I 3 (z)dz\\ 1 < C|A| 2 Yl [ dze zt (\z + ie 1 \\z + ie 2 \r 1 <Ce cx2t \\\ 2 \\og\\\\ (11.31) 

The bounds (|11.30|) and (jll.311) yield the claim. 



□ 

For completeness, we summarize Lemmata 111.11 and 1 1 1 . 21 in a single statement that is widely known as Davies' 
weak coupling limit. 

Proposition 11.1. For any t < oo, and A small enough, 

sup \\Z t - eC-^+^lao* < Ce ct |A| 2Q (11.32) 

t<A- 2 t 

11.2 Analysis of the semigroup e*® 

We will now analyze the Markov semigroup e '( _lLs + A " M ) that was derived in the previous section. The splitting of 
the generator into Ls an d A 2 M is logical from the point of view of the derivation from the microscopic dynamics, 
since M represents the contribution due to the interaction with the environment, but it is not optimal for the 
analysis of the semigroup, instead, we write 

- L s + \ 2 M = -L spin + A 2 Q, withQ:=-ad( — A) + M (11.33) 

such that on the RHS, the first term is of 0(1) and the second is proportional to A 2 . Moreover, both terms commute 
and consequently it suffices to investigate the semigroup 1 1-> e tA ® , or simply 1 1-> e t( ^ . 



11.2.1 Fourier transform revisited 

Note that up to this point, we have viewed p as a function p(x, s) on Ao and M and Q as kernels K[x' , s'; x, s) on 
Ao x Ao. In particular, translation invariant kernels were studied in terms of the Fourier transform K(p;s,s') in 
the variable x' — x. For the present discussion, it is easier to use a slightly different, in fact more natural, Fourier 
representation. Namely we replace the previously used position variable y XL + XR - J by XL + XR , Recall from Section 
12.4.11 that p can be also be viewed as a kernel p(xl,xr) £ /&(&), with Xl,Xr £ Z d . Similarly, any K in £$ can be 
viewed as a kernel 

K{x' Ll x' R , x L , x R ) £ @{@{y)). (11.34) 

Define for p e @ 2 (Jf%) 

p(p,k)= J2 ^ p ^^e ik ^- x ^p(x L ,x R ) (11.35) 

X L ,X R 

which is in L 2 (T d x T d , 3Si{S^)). Denote by p(p) := p(p, •). 

In Section f5. 3.21 we defined the Fourier transform p(p) :— p(p, •) wrt. to the variable x — [ Xl + Xr j. The two 
transforms are closely related, namely a short calculation gives 

p(p) = l p p(p) 

where l p : l 2 (S) -> L 2 (T d , 3§ 2 (y)) is the unitary map 

(l p f)(k) =^/( U ,r ? )e 1 ^e i ^ + "). (11.36) 

v . r\ 

where on the RHS, f(v,rj) £ 3§2{<5 P ), and we recall that v,rj are the (position-like) coordinates that, together with 
cl, e.R £ <j(H ap i n ), constitute the variable s £ S. 

For a translation invariant kernel K with ||if || 7 < oo for some 7 > 0, we have as before 

(Kp)(p) = K(p)p(p). (11.37) 

where K (p) £ 3§(L 2 (T d , ^(^))) (cfr. the Remark in Section l5.2.2[) . Since Q is translation-invariant and [Q, L S pm] — 
0, as we see from (|11.5[) . we can decompose further 

Q(p) = ®eesQe{p) (11-38) 

where Q e (p) ■— P E Q{p)P E with P £ the spectral projections of L sp [ n . By Assumption ICl the spaces P £ {SS{,S^)),e 7^ 
are one-dimensional and and hence we can identify 

Q 6 ( P )£<%(L 2 (T d )), e^Q 

The space Pq{SS(,S^)) has dim ^-dimensions and there is a natural basis labelled by e £ <j(H sp i n ), such that 
P (L 2 (T d , &2{y))) is identified with L 2 (T), with T = a(H spin ) x T d and so 

Qo(p)£3§(L 2 (T)). 

11.2.2 Lindblad representation for Q 

We can write a more explicit expression for M of (jll.51) using (|11.2j) and (|10.25l) : 



/ />oo />0 

(Mp)(x,y) = Y, Ue(x - y)W: P (x,y)W e ~ J dsC(0, s)<T iss W e W: p{x, y) - J 



sC(0, s)e-^W e W; P (x, y) - / d<(0, s)e~ lse p(x, y)W c W* e 

(11.39) 
In (|11.39[) we have denoted W £ :— P e (W), ( £ is the Fourier transform of ( defined in (|2.22p and we use the notation 
p(x,y) £ SS{.y) with x,y £ Z d (see Section E2|- To derive (jll.391) from (|10.25p and (jTOj) , one needs to calculate 

P £ T((W (g> l x ) a <g> U° ®{W® l y ) b )P £ = P £ W a U°W b P £ <g> l£l£. 
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To proceed write P e in terms of the projectors 7r e £ 3§(^) to the basis vector labelled by e € <j(H sp i n ). For 
7^ e = ei — e2 we have P ey o = ir ei pir e2 and for e = we have Pop = J2 e ^ep^e- Some algebra then gives the 
representation (jll.39|) . 

Denote the operator acting on p in the first term of (jlf .391) by $: 

($p)(x,v) :=J2Ux-y)w:p&yWe. (11.40) 

One can easily check that $ is a completely positive map. 

In the second term of (111. 39ft , using ("(0, s) = £(0, —s) we have 

d S e- ise C(0,s) = (l/2)C(0)+i£ £ (0) 



with t e (0) real, and in the third term one has J_ dse lse £(0,s) = (l/2)£ c (0) — ii c (0). Thus altogether we may 
write (jll.391) as 

Mp = *(p) - i(**(l)p + p**(l)) + i^Lamb, P] (H.41) 

where $* is the adjoint of <!> w.r.t. to the trace, Trp$*(0) = TrO$(p) (whereas * in W* denotes the Hermitian 
adjoint in &{&))) i.e. 

$*(l) = ^Ce(0)WW* (11.42) 

e 

and the "Lamb shift" to the Hamiltonian by 

ffLamb := VJt e (0)W e W e *. (11.43) 

In terms of the Fourier transform variables of Section 111.2.11 we have 

($p)(p, k) = J2 I Udq)W e p(p, k - q)W; (11.44) 

where the positive measure Ce(dg) is the Fourier transform of the function (, e (x), introduced in (J2.24I) . Note that 
$(p) is independent of p and hence we will write 

$ = $(p). 

Let us next decompose as in (|11.38l) . We get 

$ e = if e^O. 

Indeed, writing e = e\ — e-i, we have P e (W*,{P e p)W e i) — n ei W*,n ei pTr e2 W e nr e2 which obviously unless e' = 0. 
However, for e' = 0, £ e / = by Assumption [Cl 

Consider then <&o. As explained in the previous section it can be viewed as acting on functions <p(e,k) on 
c(-ffspm) x T d . Recalling the definition (|2.25[) of the jump rates j, eq. (|11.44[) becomes 

($ <p)(e',k , )= ]T [ 3(e',dk;e,0Me,k'-k) (11.45) 

and since Q is a finite positive measure, $>o defines a bounded positivity preserving operator on L 1 (J r ), L°°(J r ), 
and, by Riesz-Thorin interpolation, also on L q {J-), 1 < q < 00. 
Denote the escape rates of the Markov process by 

w(e):=Y,[j(e',M';e,0) (11.46) 

e' 

Then we end up with the following explicit decomposition fjll.381) : 

60 



Proposition 11.2. (a) For e — 

Qo(p) = ®o + qo(p) (11-47) 

where $o is a compact operator on L 9 (J r ),q > 1. It is also positivity improving, i.e. if ip > 0, then e'*°^9 > for 
any t > 0. qa(p) is a multiplication operator on L 2 (JP) by the function 

q ( P )(e, k) = -w(e) + iE kin (p, k) (11.48) 

where E kin (p, k) = £- 2 J= i(cos(^ + k j) ~ cos(^- - kj)) and 

mint/; := minw(e) > 0. (11.49) 

e 

(b) For e — e — e' ^ 0, Q e {p) is a multiplication operator on L q (J-) by the function 

q £ {p)(e, k) = -l/2(iu(e) + w{e')) + \E kin {p, k) + ig(e) (11.50) 

with q(e) real. 

Proof. Compactness of $o follows from the fact the function x — ¥ £ e (x) decays at inhnity for any e £ £ by Assumption 
[U] e'*° is positivity improving since by the same assumption the Markov process is irreducible. For the same reason 
the escape rates are strictly positive, i.e. (|11.49|) . 

In (|11.48[) the first term on the RHS comes from the $*(1) terms in (jll.411) and the second one from the 
Laplacian term in ([11.331) . In (lll.50[) the last term comes from the last term in (I11.41J) . This term is zero for e = 0. 

□ 

11.2.3 Spectral properties of Q £ (p) 

We can now state the main result on spectral properties of Q £ (p): 

Proposition 11.3. There are strictly positive constants 7Q,pQ,aQ,&Q s.t. the following holds. 

(a) Q e (p) extend to an analytic family of bounded operators on L 1 ^) for e = and on L l (T ) for e / for 
\lmp\ < 7q. 

(b) Take p with |Rep| < pQ and |Imp| < jq. Then the spectrum of Qo{p) lies in the half plane Rez < — oq except 
for a simple isolated eigenvalue at /q(p) with Re/g(p) > — aq/2. Moreover /q(p) is analytic and satisfies 

l/Q(p) + D QP 2 \ < C\p\ 4 (11.51) 

for a strictly positive 'diffusion constant' Dq > 0. The corresponding left eigenvector can be chosen to be analytic in 
p such that, at p = 0, it equals Ijr and the right one /J,Q(k, e) such that, at p — 0, /iq(fc, e) is constant in k, positive 
and bounded from below. The spectrum of Q £ {p) lies in the half plane Rez < — oq. 

(c) Take p with |Rep| > Pq and |Imp| < 7q. Then o~{Q e (p) C {z : Rez < — &q} for all e. 

(d) The claims (a)- (c) remain true for the operators e K k Q 6 (p)e~ K k for kgC, |k| < k for some k > 0. 

Proof. The analyticity follows since by Proposition 111.21 the only p-dependence of Q e (p)) is in the multiplication 
operators which are analytic. The spectral claim for e ^ follows from Proposition 111.21 (b) since the real part of 
the multiplication operator is positive. 

Let now e = 0. Since $o is compact and the spectrum of qo lies in a half space {z : Rez < —b} with b > 
the spectrum of Qo(p) in {z : Rez > —6/2} consists of a finite number of eigenvalues of finite multiplicity. Since 
Qo(0) generates a Markov semigroup the spectrum of Qq(0) lies in {z : Re z < 0}. Since the Markov process is 
irreducible, is a simple eigenvalue and all other eigenvalues have negative real parts. By analyticity for 7q,Pq 
small enough the simple eigenvalue persists and lies in Rez > —oq/2 whereas the rest of spectrum is in Rez < —oq 
for some oq > 0. Since the semigroup is positivity improving by the Perron- Frobenius theorem the corresponding 
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eigenfunction p,Q(k,e) can be chosen to be strictly positive. It is the density of the unique invariant state of the 
Markov process on T '. From the fact that the jump rate depends on k, k' only through k' — k, it follows that 
(iQ^k + q, e) is also invariant i.e. that fiQ is constant in k. 

Let's turn to the claim (lll.51|) . By perturbation theory of isolated eigenvalues, we get (recall that (•, •) is the 
pairing between L 1 and L°° .) 

f Q (p) = (l,(Qo(p)-Qo(0))n Q )+0(\ P \ 2 ) (11.52) 

The first term on the RHS vanishes by the explicit expression (|11.48j) and the fact that \iq is constant in k. To 
determine the second order contribution to /q(p) we invoke second order spectral perturbation theory, yielding 

d 
Iq(p) = " E ftPi<MiQo(OrW«3> +0(\p\% (11.53) 

where Vi = —id Pi E]^ n (p, •) is the "group velocity" (J2.39J) . The inverse of Qo(0) in this formula is well-defined because 
R Q (0)vifiQ = where 

R Q (p) = \HQ(p))(l\ (11.54) 

is the spectral projection to fiQ. 

The non- negativity of Dq is deduced from the fact that the operators Q p (for Imp = 0) generate contractive 
semigroups. To establish the strict positivity of Dq, we employ a standard construction: First, consider the space 
L 2 (J r , /Uq), defined by the scalar product (ip, ip'}fj, Q ■= (fiQip, t/*') where ipfiQ is the pointwise product of two functions, 
and let the operator K be defined by 

(V,^) mo ={mqOo(o)V,^) (11-55) 

where <3o(0)* is the adjoint of Qo(0) acting on L°°. By using that the density pq is bounded below and the explicit 
expression for the jump rates j (•',•), we check that the operator K is bounded and sectorial, that is, there is a 
constant Ck > such that 

|(V,(ImK)^ Q | < -C K (^,(Reif)^) w (11.56) 

where 2ReK — K + K*,21m.K = K — K* and the adjoint is taken in the Hilbert space L 2 (F,hq). Both the left 
and right eigenvectors of K are 1 (constant function on J-) and the diffusion constant can be represented as 

D Q S i j = (v i ,K- 1 Vj) MQ (11.57) 

Moreover, K inherits the unicity of the eigenvalue at and the spectral gap from Qo(0) and the sectoriality allows 
to conclude that Dq ^ 0. In the case where (3\ =02, the operator K is self-adjoint and the above reasoning can be 
slightly simplified. Hence the claims for Qo(0) hold. 

Consider then Qo(p) for p ^ 0. From (|11.48|) we see that Qo(p) is the sum of two bounded generators of 
contractive semigroups on L l (F). Hence it generates a contractive semigroup as well, as follows by the Trotter 
product formula. We need to exclude point spectrum on the imaginary axis. Thus, assume that Qo(p) has an 
eigenvalue on the imaginary axis, X p . Since the spectrum of Qo(p) is discrete in the neighborhood of the imaginary 
axis, X p must be an eigenvalue of the adjoint Q*(p), acting on L°°{F). The adjoint thus has at least one eigenvector 
ip E L°°(F). Such an eigenvector can be extended to an eigenvector of the adjoint of Q, Q*, acting on 33{Jtf?§). 
Indeed, let 

S v<p = O v e ipX G m{,^k) (11.58) 

where O^ £ £%(J#s) is the operator that acts by multiplication with <j>{e, k) on functions in j£g ~ L 2 (T) and X is 
the operator acting by multiplication with x £ Z d on Z 2 (Z d ). By definition of the Fourier representation 

Tt[^ x O <t>P }= J2 [dk<f>(e,k)(P p)(P,k,e) = (&Pop(p)) (11.59) 

e£a(H spin ) J 

where we wrote (•, •) for the pairing between functions in L 00 ^) and L 1 (J r ). Hence, for all p £ £S\{jffi$) 

Tr[S v , P Qp] = (<p, QoPoKp)) - (Q*o(P)<P> P 0P(P)) = Kin P oP(p)) = \ ^[S^pp] (11.60) 
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and it follows that Q*S V . P = XpS^^p. In Section 1 1 1 . 2 .41 we show that this equation has a non-trivial solution only 
when p = 0. 

Finally, for (c) note that the operator dk commutes with $o so that 

QZ(p) := e Kd «Q £ (p)e~ Kdk = Qe(p) + K^Unip, k - k) - E kin (p, k)) (11.61) 

Hence by perturbation theory (a) and (b) remain valid if \k\ is small enough. □ 

11.2.4 Perron-Probenius theorem for maps on BB{M£) 

Definition 11.1. A positive map <fi on £§(3%$) is ergodic if for any A > 0, A ^ 0, e^A > for any t > 0. 

Proposition 11.4 (Schrader [23] ). Assume that a positive map 4> on £§(J?s) is ergodic and that its spectral radius 
r(cj>) is an eigenvalue with corresponding eigenvector S. Then r(<fi) is a simple eigenvalue and S can be chosen 
positive definite: S > 0. 

We will apply the above theorem to the map e"^ , with t > arbitrary. First we establish 

Lemma 11.3. The map $* is ergodic on 33(,3%§). 

Proof. By duality, it suffices to prove that e**(p) > for any nonnegative p £ 8$\. Employing the fiber decom- 
position, we see that this is equivalent to e** e=0 </? > for any nonnegative function <p on T. This was proven in 
Proposition [TO □ 

Lemma 11.4. The map e t( ^ is ergodic on &(J#s). 

Proof. First, we note that the generator Q* has the form 

Q*{0)=^{0)+BO + OB\ B = i(— A + H lamb )-l$*{l) (11.62) 

Since V t : O — > e tB Oe tB is a completely positive map, for any operator B, all integrands in the Duhamel expansion 

oo „ 

° tQt = 1+ zC / dt m ---dtiQt{tt,...,t m ) (11-63) 

m=1 0<ii<...<t m <t 

with Qt(ti, . . . , t m ) := Vt-t m <&* ■ ■ ■ ^*V t2 -t 1 ^*V tl are completely positive maps, as well. 

Therefore, for any O € £$(J%), it suffices to find an too and a subset of positive measure of the TOo-dimensional 
simplex such that Qt(ti, ■ ■ ■ ,t m )(S) > 0. Let Too be smallest natural number to such that ($*) m O > 0. The 
ergodicity of $* ensures that too < oo. Then by continuity, there is an neighborhood of ti = ti = . . . = t ma = in 
the simplex where the integrand is bounded away from 0. □ 

11.3 Proof of Proposition 16.11 for n = 

Recall that To = Z x -2 io . From Proposition lll.il the relation Ls + \ 2 M — L sp i n + X 2 Q and the fact that L sp i n and 
Q, we infer 

\\T - e- u "" <oi «-»e io<3 ||20T D < C \M 2a - (11-64) 

Then we set 

T Q := c - iA_2i oi B pin e toQ = ^e" iloA_2e P £ e toQ P e . (11.65) 

£ 

We need to translate the spectral information of Proposition 1 1 1 .31 to the 7-norm. Recall that e io< ^^ = 2r 1 e to< ^ p - ) 2„ 
and that from ||11.36|l we have 
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Thus 

e^ p \v',r]';v,r]) = ^^'^-'^(I'-V™:^)^' ,rj';v,rj) (11.66) 

where we use the notation (jll.611) and both sides are operators on £$2 {•&) ■ Since 

(l; 1 e w ~*( r '2 p )(t;',t)>,f)) = ( (^,e t0, ^ p V) (11-67) 

where </?(&) = e lp ^ e lk ^ 2v+ ^ (and y/ similarly) we conclude 

\\P e e io ^ p) P e y <C sup ||e to< 5^|| (11.68) 

N<7o 

where the norm on the RHS is the operator norm on £i§(L 1 (J£ d )) 1 ^ > (L 1 (J r )), depending on e. By similar reasoning, 
we get a bound on the || • || 7 of the projection Rq(p). 

We now finally choose 70 := min(Ko,7g) where the latter constants were introduced in Proposition II 1 .31 That 
proposition can then be recast as follows. 

(1) Let p with |Rep| < pQ and |Imp| < 70. Then 

T Q (p) = e^^R Q {p) + (1 - Rq(p))T q (p). (11.69) 

with 

||(1 - R Q {p))T Q {p)\y < Ce-^\ (11.70) 

and /q(p) and Rq(p) are as in Proposition 1 1 1 .31 with 

| e t/o(p)| > ce'^-^. (11.71) 

(2) In the region |Rep| > pQ and |Imp| < 70 

\\T Q ( P )\y < Ce-^ ia . (11.72) 



We can now apply spectral perturbation theory, e.g. Lemma lA. II) . First, for to sufficiently large (compared to the 
constants in the above statements), the eigenvalue e to ^ Q ^ is isolated, by (|11.70[) and (|11.71l) . Thus, taking then A 
small enough, |A| < A(to) the following holds. 

(a) The isolated eigenvalue persists and the new eigenvalue e^ p ' has the same properties: /o(0) = by unitarity, 
V/o(0) = and \fo{p) + D p 2 \ < C\p\ 4 by lattice symmetries, with Dq -> toDQ as A — >• 0. 

(b) The bounds (|11.70[) and (|11.72p hold for e /o(p) and T (and different C and c). Now take r large enough and 
then Proposition 16. II (except for (J6.12I) and (I6.4J) for T n (0,p)) holds for n = and to £ [to,2to]. 

To get the claims for Ro{0), i.e. (16.121) . note that (|10.53[) tells us that /x' 3 is A 2 -close to the A = system 
equilibrium state, i.e. the Gibbs state ~ e _/3i/Bpin . The latter is identical to the state £*q(0) in case fi\ = $% = (3, 
and hence one has \\PP - Rq(0)\\ w = C(A 2 ). Since \\Rq(0) - Rq(0)\\ = C(|A| 2a ), (l6~12l holds for n = 0. Finally, 
the bound (|6.4[) for T o (0,p) was proven in Section [10.31 

A Appendix: Spectral perturbation theory 

In the lemma below, we gather some spectral perturbation theory that is used in the proof of Lemma 17.31 We 
consider operators on a Banach space &/. Denote by || • || the usual operator norm on 38{stf) and let || • || by a norm 
on (a dense subset of) £$(&/) such that 

Mil < ML 

\\AB\l < \\A\U\B\l 

The following lemma could just as well (and more naturally) be stated for || • || as for || • || . The use of || ■ || is 
dictated by our application, where srf = 1°°{S) and || • || = || ■ || 7 . In fact, in Lemma T7.3I we did not mention the 
space srf — l°°(S) at all, but note that without such an underlying space, it is not a-priori clear what one means by 
expressions like 'a simple eigenvalue'. 
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Lemma A.l. Let the operator A$ have a simple eigenvalue ao with corresponding one- dimensional spectral projection 
Po such that 

\\A - aoPollo < |a |, ||P || <> <oo (Al) 

We study A = Aq + A\ for some perturbation A\ . Define 

6(r):= ^k + 1 /or r<|a |- || Ao-aoPoL (A2) 

r |ao|-r- H^o-aoPollo 

//, for some r > 0, we have 

||i4i|| 6(r) < 1 (A3) 

then the eigenvalue persists when the perturbation is added (we call it a) and 

\a-a \<r (A4) 

The spectral projection P corresponding to a satisfies 

ii P -P | k < 2Trt(f) J^Mr)_ (A5) 

Proof. The condition (jAlj) implies that the eigenvalue ao is an isolated point in the spectrum. Let C r C C be the 
circle with radius r centered on ao and such that Aq has no other spectrum within C r (the latter condition will be 
implied in our case by ||j4o ~~ to Poll < r). To prove that the eigenvalue persists within the circle C r , it suffices to 
show that sup 2SC \\(z— (Ao + kAi))~ 1 \\ remains finite for k € [0,1]. Without loss, we can assume k = 1 in the proof 
below. Consider the Neumann series 

(z - A)- 1 = (z- Ao)" 1 J2 (M* - AarT (A6) 

n>0 

To establish that the resolvent remains finite, we show a stronger property (because of the different norms), namely 

sup^(||A 1 |U||( z -A )- 1 |Ur<oo (A7) 



To get this, we bound 



,-I.M . .... ||PO||0 , . , • ... , ,,,., 1 



sup \\(z - A )~ l \l < sup " u "°. + su P \\(z-(A -aoPo))- l \\ <b(r) (A8) 

zec r zec r V- — a o| zee,. 



where the last equality follows from the Neumann series. Therefore, (|A7j) is finite if ||Ai|| o(r) < 1 and we conclude 
that the eigenvalue a lies within the circle C r . The bound for the projection follows as 

P-Po = / dz(z - A)- 1 - (z - Ao)- 1 (A9) 

Jc r 

= j dz(z-^o)- 1 (j2A 1 (z~A a )- 1 \ (A10) 

11^-^olU < / dziKz-Ao)- 1 !! "tl'fr^^-iH (ah) 

Jc r i-ll^illoll^-^o) lo 

< Mr) ' 1 ^! ( Al2 ) 

n 
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Let us denote by c, C constants that depend only on ||Po||©- If we assume that 

IMli|L<c(|ao|-||4)-aoPoU, (A13) 

then we can choose r = C||^4i || to estimate the spectral shift. To bound the difference of projections, we choose 
r — c(\ao\ — \\Aq — aoPollo) an d we obtain 

|p - f * sc RHFMg (A14) 

B Appendix: Combinatorics 

We collect some easy combinatorial lemma's. 

Lemma B.l. Let ^ be a spanning tree on the set A c No, and define 

dist(^):= JJ (1+|t'-t|) (Bl) 

{r,T'}eS(sr) 

where £"(2?) is the set of edges of & . Then 

dist(A) < dist(^) (B2) 

with dist(A) defined as in Section^ 

See e.g. [22] for the 5-line proof. The following result lies at the heart of cluster expansions. For subsets S £ No, 
write S ~ S' whenever S S' 7^ 0. Let S be a collection of subsets S = {Si, . . . , S m }- Note that the connectedness 
of the graph with vertices Sj and edges {Si, Sj} whenever Sj ~ Sj is equivalent to the connectedness of S defined 
in Section |9] and, as in that section, we write £ for the set of connected collections. 

Proposition B.l. Let w(-) be a function on subsets S G No- Assume that for some k > 0, 

J2 e"W\w(S)\<K\S'\. (B3) 

S:S~S' 

Then, with w(S) = Yises w (^)> 

^Us^useeM^l < e KlS ' 1 (B4) 



We refer to [25] for the proof. The condition (JB3I) is often called the Kotecky-Preiss condition. An immediate 
consequence is 

J2 HS)\< £ | W (5')|£l [m u5 e c]K5)|< Y, HS')\e KlS ' 1 < k (B5) 

5e£:Supp59r S':S'3t S S':S'3t 
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